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Transmission of Sound through Air at Low Pressure 


R. B. Linpsay 
Brown University, Providence 12, Rhode Island 


HE famous elementary physics demonstra- 

tion of the electric bell ringing in a bell 
jar from which the air is gradually pumped out, 
and the corresponding decrease in the intensity 
of the sound that is heard, is familiar to all. 
We tell our students that this experiment proves 
that a material medium is necessary for the trans- 
mission of sound—one of those glib statements 
which make one wonder whether it is ever 
possible to tell the truth about anything in an 


elementary physics course. Actually the experi-_ 


ment demonstrates the impedance mismatch be- 
tween transducer and surrounding medium, but 
an understanding of this explanation demands a 
knowledge of the theory of wave propagation. 

This paper is limited to a survey of the at- 
tenuation of sound in air as a function of pres- 
sure or density. For the purpose of unification 
and simplification we shall treat the problem in 
terms of relaxation dissipation. 


Wave Propagation with Relaxation Dissipation 


The usual treatment of relaxation dissipation 
in compressional-wave propagation leads to a 
complex wave velocity which, when substituted 
into the ordinary harmonic wave solution of the 
standard wave equation, provides for the ap- 
pearance of the dissipative term. It seems worth 
while to present here a somewhat different point 
of view, which removes some of the mystery 
associated with the complex velocity and also 
emphasizes the fact that the equation we have 


to deal with in the description of attenuation is 
not the standard wave equation anyway. 

For simplicity we shall confine our discussion 
to plane compressional waves in one dimension. 
If we denote the displacement of the elastic 
medium in the x direction by £, the excess pres- 
sure by p, and the mean density by po, the usual 
hydrodynamic equation of motion for small de- 
partures from equilibrium is 


po(0°E/ A) = — (dp-/9x). (1) 


In the common method of deriving the wave 
equation, Eq. (1) is supplemented by Hooke’s 
law in the form 


be = Bp./po= — B(dé/dx), (2) 


where p, is the excess density and B the bulk 
modulus. Elimination of p, between Eqs. (1) 
and (2) yields the well-known equation 


at/ar=c(a%/ax), (3) 
where the wave velocity is 
c=(B/po)}. (4) 


Now Eq. (2) is a static-equilibrium equation 
which ignores the dynamical mechanism by 
which the density is changed through the agency 
of the excess pressure. To take account of the 
latter we must replace this equation by the more 
accurate pressure relaxation equation 


Pe=Bp/potRpe, (5) 


where partial differentiation with respect to 
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time is indicated by the use of the dot. The 
quantity R may be called the relaxation con- 
stant. The corresponding relaxation time is 
poR/B. This is the time taken by the density 
to attain to the fraction 1—1/e of the equi- 
librium value corresponding to a given imposed 
excess pressure. If we employ p.= —po(d£/dx), 
we can rewrite Eq. (5) in the form 


















(6) 


Elimination of p, between Eqs. (1) and (6) now 
leads to the equation 


Oe ort arg 
(B / Po) — +R— 


Ox? ax? 


= (7) 
ar 


which is formally equivalent to the Stokes equa- 
tion for the transmission of a plane compres- 
sional wave through a viscous medium! if R is 
treated as proportional to the viscosity coeff- 
cient pu. Indeed, in the standard treatment of 
viscosity dissipation, 


R=4y/3po. (8) 


However, for the moment let us proceed 
without too great specialization. It will be 
sufficient to consider a harmonic particular solu- 
tion of Eq. (7) in the form 


§=exp(—ax) cos(wt—kx). 


(9) 


1See Lord Rayleigh, Theory of sound, (Dover Publica- 
tions, 1945), sec. 346, Eq. (4). 
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Fic. 1. Classical approxima- 
tion: av<1; 2a(db/cm) as a 
function of frequency for various 
pressures. 


Substitution shows that this is a solution of 
Eq. (7) provided 


wT po 
ak =——_—_—__, (10) 
B(i+w?r?) 
and 
w*po 
k? — o? =—____—__, (11) 
B(1+w?r?) 
where 
T=poR/B (12) 


is the relaxation time. The elimination of k 
between Eqs. (10) and (11) yields for a, which 
is the amplitude-attenuation coefficient, the bi- 
quadratic equation 
wpo w®72 p¢2 
a4 -+——_—_____,?—--—____—_——-=Q._ (13) 
B(1+w?r?) 4B?(1+w?*r?)? 


The real solution is 


1 4 6 (1+w?7?)!—17 
2 PN ig 
2B/po | (1+ w?7?) 
This is in rather awkward form. If we solve for 
k from Eq. (10) we get 
‘ wr { (1+ wr?) [ (1 +7”)! — 1 ]}-? 
(2B/po) 


But this is w/V, where V is the phase velocity 
of the wave. Hence 


(2B/po)* 
GF ghiinens 


QT 





(15) 


((1+0%7)[(1+e%)'—-1}}4. (16) 
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lf now we express a in terms of V, the result is 
wr V 
(2B/po)(1+-u%r2) 


Since V depends on w—that is, there is dis- 
persion—Eq. (17) does not exhibit directly the 
dependence of a on the frequency. The ultimate 


approximate form will depend on the magnitude 
of 7. 


(17) 


Viscosity and Heat-Conducticen Dissipation 


Viscosity dissipation can be treated by the 
method of the preceding section. Here, from 
Eqs. (8) and (12), 


7=4y/3B. (18) 


By simple kinetic theory, u» is independent of 
the pressure. On the other hand, the adiabatic 
bulk modulus is yp, where p is the pressure and 
y the ratio of the specific heat at constant pres- 
sure to that at constant volume. Since for air 
at room temperature yu is 1.8 10-4 dyne sec/cm?, 
it follows that for pressures above 1 dyne/cm? 
and even considerably below, the value of 7 in 
seconds is very small compared with unity. 
For not too high frequencies we therefore have 
wri. Under these circumstances Eq. (17) re- 
duces to the conventional form 


a = 2yw*/3poc’, (19) 


where 
c=(B/po)}. (20) 


The effect of heat conduction also is, to a 
first approximation, a relaxation effect. We shall 
not discuss it here. It is well known that it 
leads, for sufficiently low values of thermal con- 
ductivity and not too high frequencies, to an 
absorption coefficient a of the form 


w (—): 
2poce\ » Jey 
where 7 is the thermal conductivity and c, the 


specific heat at constant volume. The corre- 
sponding value of the relaxation time becomes 


y-—1 
-(— bee 
Y 


2 See Lord Rayleigh, reference 1, sec. 349. 


(21) 


(22) 
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It is of interest to observe that for air at normal 
atmospheric pressure and room temperature, the 
value of 7 calculated from either Eq. (18) or 
Eq. (22) is of the order of 10~° sec. 

We now introduce the standard kinetic theory 
results, 


(23) 
(24) 


where / is the mean free path and »,, the root- 
mean-square molecular velocity. Using the for- 
mulas im=(3p/p)? and c=(yp/p)!, the latter 
being sufficiently accurate for an approximate 
calculation, we get for the intensity absorption 
coefficient, taking into account both viscosity 
and heat conduction, 


n= 2com, 


p=0.3 porn, 


_ 23.6041 
(25) 


(3/y) (-¥ i 


If we specialize at once to air with y=1.41 and 
c=3.44X10* cm/sec at 20°C, this becomes in 
reciprocal centimeters 


2a=2.77v71X 10-8, (26) 


if v is frequency in cycles per second and 7 is 
the mean free path in centimeters. For the mean 
free path we use the standard Clausius formula, 


l=3m/4rD" po, (27) 


where m is the mass and D the diameter of the 
molecule. Using D=3.1X10-* cm as the diam- 
eter of the nitrogen molecule (that of oxygen is 
about 3X10-* cm) and m=46.5X10-" g, we 
obtain for the mean free path, if po is in grams 


Fic. 2. Classical approximation: a\ <1; 2a (db/cm) as a 
function of frequency for various pressures. 
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per cubic centimeter, 
l= 1.15 X10-8/po, (28) 

and 2a becomes, in reciprocal centimeters, 

2a = 3.19 X 10-16 y?/ pp. (29) 


In more conventional units the transmission loss 
in decibels per centimeter is 


Qa = 13.9 10-"%r2/po. (30) 


It must be recalled from (17) coupled with wr<1 
that the limit on the validity of Eq. (30) is 
given by a\<1 or 


5.5 10-"v/po <1. (31) 


In Fig. 1 we present a logarithmic graph of Eq. 
(30) for various air pressures, from 760 to 107! 
mm of mercury. Figure 2 extends this to pres- 
sures from 10~? to 10~* mm of mercury. The fre- 
quency limits to the validity are indicated by 
the length of the line in each case. Thus at 
atmospheric pressure, Eq. (30) is approximately 
valid for »<10° c/sec. Obviously the upper 
frequency limit for validity decreases in the same 
proportion as the pressure. Except for sound of 
frequency less than 1 c/sec, Eq. (30) has no 
validity for pressures lower than 10-5 mm of 
mercury. 

It is interesting to note that the experimental 
values of absorption in dry air at atmospheric 
pressure and 26.5°C quoted by Sivian* corre- 
spond to 


2a = 19.7 K10-'*y?/ po (db/cm), 


with a good frequency-squared variation from 
20 to 200 kc/sec. The increase over the theoreti- 
cal value given by Eq. (30) is not one order of 
magnitude and is not sufficient to show up on 
the scale used in Fig. 1. 


Calculations of Tsien and Schamberg 


Recently Tsien and Schamberg* have em- 
ployed the second-order approximation to the 


3L. J. Sivian, J. Acous. Soc. Am. 19, 914 (1947). 

4H.S. Tsien and R. Schamberg, J. Acous. Soc. Am. 18, 
334 (1946). As pointed out by C. S. Wang Chang and G. 
E. Uhlenbeck, there exists an error in the numerical values 
taken by Tsien and Schamberg from Chapman and Cowl- 
ing (reference 5). This led to incorrect values for bs and 
bs in Table I of Tsien and Schamberg’s article. When cor- 
rected these become b;= —2.29, b;=+9.05, and these 


values are used in the present review to calculate 2@ in 
Eq. (32). 
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viscous stresses and the heat flow given in the 
more recent kinetic theory considerations of 
Chapman and Cowling’ to calculate the absorp- 
tion of plane sound waves in rarefied mediums. 
They obtain an expression for 2a in the form of 
a power series in v?. When the values for dry air 
at 20°C are inserted the series becomes 


2a =9.04 X 10-1(v2/p)(1—9.15 K 10 y2/ p? 
+1.34X10-p!/pi+---). (32) 


Here # is the pressure in millimeters of mercury. 
Reference 4 shows that the convergence is as- 
sured only if the Reynolds number 


R’ =1.64X10°p/v>1, (33) 
or 


p/p <1.64X 10°. 
In terms of density po, this condition becomes 
0.97 K 10-" v/ po <1. (34) 


For given po this provides a somewhat greater 
range of frequency validity than the simple 
classical theory as exemplified in Eq. (31), 
though the order of magnitude is the same. As 
a matter of fact, the contribution of the correc- 
tion terms in Eq. (32) in the region of validity 
is very small, too small indeed to show up in 
Figs. 1 and 2. 


Classical Approximation for Very 
High Frequency 


The relatively small magnitude of the correc- 
tion introduced by the considerations of Tsien 
and Schamberg encourages us to examine once 
again the classical expressions for @ for ultra 
high frequencies at ordinary pressures or ordi- 
nary frequencies at very low pressures. It is 
indeed too much to expect that these expressions 
will have more than order-of-magnitude sig- 
nificance, and further work along the lines of 
the theory of nonuniform gases is definitely 
indicated. If we revert to Eq. (14) we see that 
for wr>1, a assumes the very simple form 


a=2rv(2we?r)-}, (35) 
where for air, using the data already given 


wr = 1.23 10-"»/ po. (36) 





5S. Chapman and T. G. Cowling, The mathematical 
theory of non-uniform gases (Cambridge Univ. Press, 1939). 
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Fic. 3. Comparison between 
low and high frequency approxi- 
mation for the absorption co- 
efficient. Frequency range, 10% 
10° c/sec. 


For v>107 c/sec and p<1 mm of mercury the 
conditions leading to Eq. (35) are met sufficiently 
well, and the same is true for v>10? c/sec and 
p=10-5 mm of mercury. Under these conditions 
a varies effectively as v?. 

This treatment, of course, neglects the effect 
of heat conduction. The simplest way to intro- 
duce this effect in the general case is to go back 
to the treatment by Rayleigh,® where it is shown 
that owing to the combined effects of viscosity 
and heat conduction the particle velocity is 
given by 


u~exp(—x(2:)4), (37) 


where 2 is the smaller of the roots of the 
quadratic 


2°(c?/y+4iwpu/3po) 
— 2? +10(n/Copo+4u/3po) ] 


X te vpo/n —1w*Copo/n=0. (38) 


The various symbols in this equation have the 
same meaning as earlier in the paper. If we solve 
the quadratic under the approximation 


pw/ po>c?, (39) 


which is valid at high frequencies and low pres- 
sures, we obtain as a sufficient approximation 
to the required root 


21 = (twpo/2u)(1+0.42poc?/ pw). (40) 
On this basis the absorption coefficient becomes 
2a = 8.04 X 10?(pov)?. (41) 


to an approximation that neglects poc?/yw in 
comparison with unity. 


6 Reference 1, secs. 348 and 349. 


The values of 2a from Eq. (41) are plotted in 
their appropriate regions of validity for various 
pressures in Fig. 3. For purposes of comparison, 
the classical low-frequency approximation values 
(from Figs. 1 and 2) are also plotted on the same 
graph. No attempt has been made to interpolate 
between the two regions. Of interest is the fact 
that at the higher frequencies the absorption 
coefficient varies as the square root of the fre- 
quency and also as the square root of the density. 
Hence the curves appear in reverse order as the 
pressure is decreased. 

The results just obtained must not be taken 
too seriously as the following considerations 
make clear. Indeed, Eq. (41) implies that, for 
given v, a0 as po—0, a manifestly absurd 
result. We certainly must not employ Eq. (41) 
for too low densities. The difficulty encountered 
here results from the assumption that the vis- 
cosity is independent of the pressure no matter 
how small the latter is. This is probably true for 
pressures as low as 10-? mm of mercury, but 
hardly for much lower pressures, as a reference 
to any standard text on kinetic theory of gases 
will confirm. As the pressure is decreased, the 
number of molecules per unit volume ultimately 
becomes too small to provide the transfer of 
momentum responsible for viscosity, and hence 
the viscosity itself will have to decrease. Pre- 
sumably under these conditions y» will fall off 
at least as fast as the first power of po. The effect 
of this would be to keep 7 relatively constant as 
po decreases and hence to remove the po from 
Eq. (41), 2@ still varying approximately as v}. 

However, let us contemplate pressures for 
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Fic. 4. (a) Schematic diagram of apparatus; (b) Orienta- 
tion of microphone and loud-speaker. 


which y is still practically constant and fre- 
quencies high enough that wr>1. Considering 
viscosity only for simplicity (these are order-of- 
magnitude considerations), we readily see from 
Eqs. (28) and (36) that 


wr~10-*VI/X, (42) 


where V is the phase velocity and X is the wave- 
length. From Eq. (16) we have, when wr>1, 


V =c(2wr)!. (43) 
This finally leads to 
wr = 25P/d (44) 


for air. Now when /~ X, we must conclude from 
elementary considerations, as was pointed out 
in 1917 by Schrédinger,’? that the propagation 
of sound becomes impossible, since under these 
conditions it is on the average impossible to 
maintain differences in density and pressure in 
neighboring parts of the medium in the face of 
the increased rapidity of the exchange of mole- 
cules accompanying the lengthened mean free 
path. It is clear that there is only a very small 
pressure frequency range for which /<\ and 
wr>>1. Hence the values of 2a given in the high- 
frequency region in Fig. 3 are probably unre- 
liable, except in the sense that they confirm the 
very large attenuation in this range which the 
common-sense view dictates. 


7E. Schrédinger, Physikal. Zeit. 18, 445 (1917). 
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Fic. 5. Tank experiment. Transmission loss (in db) 
versus pressure. 
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It would appear to be desirable to investigate 
the variation of viscosity with pressure at very 
low pressures in large spaces, that is, relatively 
unconfined mediums. The upper atmosphere im- 
mediately comes to mind in this connection. 
Most of the work on viscosity and heat conduc- 
tion has so far been restricted to relatively small 
vessels in which surface effects ultimately pre- 
dominate. 


Measurement of Acoustic Transmission Loss 
in a Low-Pressure Tank 


Few experimental data have as yet appeared 
in the generally available literature on the trans- 
mission of sound through air at low pressure, 
though several projects are underway and pre- 
liminary results have been reported recently.* 
We shall describe in this section the first results 
of an experimental program being initiated at 
Brown University. 

In a metal tank of cylindrical shape (75 cm 
in diameter and 50 cm high), used originally as 
a metallic evaporation chamber, a loud-speaker 
and a microphone were set up very close together 


H.C. Rothenberg and W. H. Pielemeier, Pennsylvania 
State College, ‘“‘Preliminary report on absorption in air 
from 10 to 100 it manned by pulse techniques,” and 
L. P. Delsasso, R. W. Leonard and R. Halley, University 
of California at Los Angeles, ‘ ‘The absorption of sound in 
air at reduced pressures;’’ papers delivered at the Wash- 
ington meeting of the Acoustical Society of America, 
April 1948. 
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(Figs. 4a and 4b). Care was taken to insulate 
the mounting as far as possible from the walls 
of the tank. The loud-speaker was driven by a 
variable-frequency oscillator. In these prelimi- 
nary experiments low frequencies—that is, under 
7 kc/sec—were used exclusively. The chamber 
could be evacuated by a mechanical pump to a 
pressure of about 5 mm of mercury. 

If we let J, be the measured intensity of the 
sound received at the microphone (in arbitrary 
units) for given loudspeaker velocity amplitude 
at atmospheric pressure and standard room tem- 
perature, and J, be the corresponding intensity 
for the same loudspeaker velocity amplitude at 
pressure p, we shall define the quantity 


10 logi(Ia/I»p) (45) 


as the transmission loss in decibels at pressure 
p relative to atmospheric pressure. The experi- 
mentally determined values of this quantity at 
various frequencies from 384 to 6950 c/sec are 
plotted in Fig. 5. The solid straight line repre- 
sents the quantity 


20 logio(poa/ pop), (46) 


where poa is the density of air at atmospheric 
pressure and po, the density at pressure p. The 
basis for expression (46) is simply that the radia- 
tion resistance for any transducer radiating into 
a fluid medium of density po is directly propor- 
tional® to po. Hence other things remaining the 
same, the intensity of radiation varies directly 
with the mean density and hence the mean 
pressure of the medium. By the reciprocity 
theorem the receiving microphone obeys the 
same rule, which accounts for the factor 20. If 
we like we may consider the quantity in the 
expression (46) as representing the theoretically 
expected loss in the pressure chamber due to 
the impedance mismatch between transducer 
and air. As the pressure of the air is reduced it 
becomes progressively harder to get the energy 
from the transducer into the air and likewise 
from the air to the microphone. 

From the general order-of-magnitude agree- 
ment between the experimental results and the 
curve based on expression (46) it seems clear that 
the pressure-tank measurements here reported 


°G. W. Stewart and R. B. Lindsay, Acoustics (Van 
Nostrand, 1930) p. 247. : 


indicate primarily the increasing impedance mis- 
match with decreasing pressure. The indicated 
shift in slope with frequency is suggestive but 
hardly significant at the present stage of the ex- 
periments. At the pressures and frequencies here 
considered actual absorption due to viscosity 
and heat conduction can hardly play much of a 
role, and this indeed is confirmed by the theo- 
retical values already noted, where for v= 7000 
c/sec and p=10 mm of mercury the loss comes 
out to be about 4.5 X 10-* db/cm. 

To use the tank technic for the study of genu- 
ine absorption of sound at low pressure it will 
be necessary to go either to much higher fre- 
quencies or to much lower pressures or to both. 
This involves the construction of a new sound 
source and receiver outfit and the attachment 
of a faster, high-power pump. 

The employment of the method also involves 
some theoretical points which must be examined 
with care for their effect on the results. These 
are the change in the source and receiver im- 
pedance due to the change in the air load at 
lower pressure and the effect of standing waves 
in the chamber. No account has so far been taken 
of the normal modes of the chamber. 

Many interesting problems remain to be at- 
tacked, among them the influence of water vapor 
and other gases in air at low pressure. This 
will introduce a consideration of the extension 
to low pressures of the Knudsen-Kneser thermal- 
relaxation-time theory of sound absorption. Pre- 
liminary examination indicates that the lowering 
of the pressure will increase the expected ab- 
sorption from this source. The transmission of 
sound waves of finite amplitude at reduced pres- 
sure will provide another challenging problem. 
The next few years will doubtless witness a vast 
increase in our knowledge of the propagation of 
high-frequency sound through highly rarefied 
mediums. This must go hand in hand with theo- 
retical studies based on the kinetic theory of 
gases at extremely low pressure. 

It is a pleasure to acknowledge the assistance 
of my colleagues Professor J. R. Frederick and 
Messrs. Carl E. Adams and E. F. Smiley in 
connection with the experimental measurements 
reported in this section, as well as certain of the 
calculations and curve plotting in the rest of the 
paper. . 


























































































































NUMBER of problems dealing with the 
physics of the earth’s upper atmosphere, 

and many problems in the field of astrophysics, 
involve a quantity which is loosely called ‘‘tem- 
perature.’’ A good deal of confusion exists re- 
garding these temperatures, and it is important 
to realize that the word may have different mean- 
ings in different cases. It is the purpose of this 
article to explain the significance of the various 
temperatures quoted in astrophysical problems. 
The concept of temperature in an atmosphere 
is complicated by departures from thermal equi- 
librium. If thermal equilibrium existed, the prob- 
lem would be much simplified and all methods 
of computing temperature would give the same 
result. Consider what is meant by thermal equi- 
librium. Suppose we have an assemblage of mole- 
cules, atoms, and ions at a given temperature T. 
Collisions will take place between particles, and 
if the energy transferred during collisions is 
sufficient, molecules and atoms may be excited 
to higher energy levels, and a short time later 
unload this energy in the form of radiation. 
Other molecules and atoms may absorb the 
radiation; hence, the energy of the assemblage is 
being transferred from particle to particle 
through collisions and absorption and emission 
of energy. Stated more fully, atoms (or mole- 
cules) in a certain energy level A may absorb 
radiant energy and be raised to a higher level B. 
After a short interval these atoms (or molecules) 
return to level A and emit energy. Atoms (or 
molecules) in the lower level may also reach B 
by collision with other particles. In this case, 
the colliding particles rebound with lessened 
energy. Furthermore, atoms (or molecules) in 
the higher level may return to A by collision 
with particles. These particles rebound with 
increased energy. If the assemblage is in thermal 
equilibrium, then the number of absorptions 
from A to B equals the number of emissions from 
B to A, and the number of collisional excitations 
from A to B equals the number of collisional de- 
excitations from B to A. When this condition 
exists no spectral lines are formed, since energy 
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levels are populated and depopulated at the 
same rate by both processes. Hence, the detec- 
tion of spectral lines in an atmosphere is an 
indication of departure from thermal equi- 
librium. 

There are a number of useful equations which 
may be applied when thermal equilibrium pre- 
vails. One of these is the familiar Maxwell 
velocity distribution law which gives the frac- 
tion of particles having velocities in a range 
between v and (v+dz), 


m 


2rkT 





3/2 
f(o)do=4r( ) v exp(—3mv?/kT)dv. (1) 


In certain cases we may have a volume of gas 
which initially is not in thermal equilibrium. 
Collisions between the gas particles will gradu- 
ally bring about a Maxwellian velocity distribu- 
tion at a rate that depends upon the pressure, 
temperature, and the nature of the particles. 
For example, the diffuse material in interstellar 
space at a space density of the order of 3X 107° g 
cm~’, will reach a Maxwellian distribution in 
about 4X 10° years. The same reasoning may be 
applied to the aggregation of stars known as the 
galaxy. The space density of stars in the neigh- 
borhood of the sun is 410-57 cm~-', and stars 
the size of our sun will suffer an actual collision 
once every 2X10!" years. Hence, for these stars 
a Maxwellian velocity distribution will be 
reached only after a long period of time of the 
order of 7X10!* years. 

The Boltzmann distribution law, which gives 
the distribution of atoms or molecules in various 
energy levels, may be applied when thermal 
equilibrium conditions exist. The law may be 
written 


Wn 
Na=No— exp(—xn/RT), (2) 


0 


where N, is the number of atoms or molecules 
in level 2, No the number of atoms or molecules 
in the ground level, w, the statistical weight of 
level n, wo the statistical weight of the ground 












TEMPERATURES DERIVED 


level, xn the excitation potential of level m and 
T the absolute temperature. 

In an atmosphere consisting of atoms, mole- 
cules, positive ions, and electrons, only the elec- 
trons may be expected to approach a condition 
of thermal equilibrium. The reason is that atoms, 
molecules, and positive ions may absorb incident 
radiation and have an energy level distribution 
which is non-Boltzmann in character. Free elec- 
trons may experience elastic collisions with other 
particles, inelastic collisions which result in the 
excitation of atoms or molecules, collisions of 
the second kind with atoms or molecules in 
excited metastable levels, free-free transitions 
in the neighborhood of nuclei with the resulting 
emission of energy, or capture by ions with the 
emission of energy. Bohm and Aller! have con- 
sidered the motions of the free electrons in the 
atmosphere of a planetary nebula, and conclude 
that ‘the process of elastic scattering is by far 
the most probable. If this is so, the free electrons 
will reach a velocity distribution which is close 
to a Maxwellian distribution. 

A third expression that applies to thermal 
equilibrium conditions is the Saha ionization 
equation, giving the relative numbers of ionized 
and neutral atoms (or molecules) in an as- 
semblage at temperature JT. Ionizations are 
brought about by collisions and absorptions of 
energy at a rate depending upon the temperature 
and the type of atoms (or molecules) ; recombina- 
tions take place at a rate depending upon the 
pressure. High temperature with low pressure 


favors the ionized condition. The Saha equation 
is 


Ma 2B,(T) (24m)?!(kT)5!2 
No BolT) hs 


exp(—x/kT), (3) 
where N, is the number of ionized atoms, No the 
number of neutral atoms, P the pressure of the 
free electrons, and x the ionization potential. 
The function Bo(T) gives the distribution of the 
neutral atoms among the various energy levels. 
It can be computed for any temperature. Fol- 
lowing the same notation, B,(T) is the corre- 
sponding function for the ionized atoms. 

These equations which have been discussed 
briefly are used to determine temperatures of 


1D. Bohm and L. H. Aller, Astrophys. J. 105, 131 (1947). 
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stellar atmospheres, of the earth’s atmosphere, 
and of the interstellar material. It is more satis- 
factory, of course, to determine temperatures by 
measuring the total radiation emitted by a star, 
or by studying the intensity distribution of this 
radiation. Temperatures calculated by the first 
method are called “effective temperatures,” 
since it is this quantity that is effective in pro- 
ducing the observed radiation from the star. As 
an example, we will now indicate how to compute 
T, the effective temperature of the sun. Stefan’s 
law gives the total radiation per unit time from 
a unit area of the solar surface, or, in the usual 
notation, E=oT‘, where o has the value 5.672 
X10->5 erg cm deg‘ sec. The amount of 
solar radiation falling normally on unit area of 
the earth’s surface per unit time is called the 
solar constant and has the value 1.35X10° erg 
cm~ sec~!. The total radiation emitted by the 
sun in unit time is 47R?F, where R is the distance 
from sun to earth and F is the solar constant. 
If 7 is the radius of the sun, its surface area is 
4rr*. Then by Stefan’s law, 


(49R?F/4ar?) =oT*. 


Substituting the known values of R, F, 7, and o, 
we find that T is around 5750°K. This method 
can be used only for those stars whose distances 
and sizes are known, data which are available in 
relatively few instances. In a few cases effective 
temperatures may be obtained from a study of 
the spectra of double stars. 

A study of the intensity distribution in the 
spectrum of a star yields a quantity known as 
the “color temperature.”’ This temperature is 
obtained by comparing the observed spectral 
energy curve with theoretical curves computed 
from Planck’s law. The theoretical curve repre- 
senting the equation which includes the stellar 
temperature will match the observed curve. This 
method of obtaining temperatures has limited 
applicability, because only in relatively few in- 
stances is it possible to determine the true spec- 
tral energy curve of a star. The continuous spec- 
trum of all distant stars is distorted by the red- 
dening effect of interstellar dust, and accurate 
corrections are difficult to make. A recent at- 
tempt? to determine the color temperatures of a 


2W. Petrie, Public. Dom. Astrophys. Observ. 7, No. 25 
(1948). 


















































































































































































380 


number of Wolf-Rayet stars has encountered 
this difficulty. 

In the absence of information on effective or 
color temperatures, it is sometimes possible to 
determine the temperature by applying Eqs. (2) 
or (3) to the intensities of certain emission lines 
in the spectrum of a star. For example, Eq. (2) 
may be used to compute a quantity called the 
“excitation temperature.”” The emission in ergs 


per cm® per sec in a spectral line is given by the 
expression, 


Enw _ N,A nn'DVan'; 


where JN, is the number of atoms in the upper 
level involved in the production of the spectral 
line, Ann is the Einstein transition coefficient 
between levels m and n’, and vn,’ is the frequency 
of the radiation associated with a transition 
between levels » and n’. If thermal equilibrium 
conditions exist we may substitute the Boltz- 


mann equation in the above expression. We then 
have 


Wn 
Enn' = No— exp(—Xxn/RT)Ann bran (4) 


0 


Taking the ratio of the intensities of two spectral 
lines from the same atoms in the same condition 
of ionization, we may write 


Ex, wAin 
—=— — — exp((x2—x1)/RT). (5) 


Substituting measured values of EZ, and E, and 
the quantities w, A, v, and x from atomic data, 
we may determine the temperature T. This 
quantity is called the excitation temperature, 
since it refers to the absorbed radiation or the 
colliding particles which excite the atoms to the 
various energy levels. Aller* has used this method 
to compute temperatures in the atmospheres of 
a number of Wolf-Rayet stars, and Petrie* has 
determined the excitation temperature of the 
earth’s upper atmosphere from the intensities of 
permitted oxygen lines in the auroral spectrum. 

If the intensities of a number of lines are 
available, it is better to proceed in the following 
manner. Write Eq. (4) in logarithmic form and 
let (No/wo)h=Ci. Then, logEnn» =logCit+logwn 


3L. H. Aller, Astrophys. J. 97, 135 (1943). 
4W. Petrie, Can. J. Research A25, 293 (1947). 
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— (5040xn/T)+log Ann +log yan. In this form 
of the equation x, is in electron volts. Now 
let waAnnYan' =D. Then log(Enn/D)=logC, 
—(5040x,/T). A plot of logE,,-/D against x, 
produces a straight line, the slope of which gives 
the temperature. Goldberg’ and Petrie® have 
used this method to determine the temperature 
of the solar chromosphere. However, if the 
Einstein A values for the spectral lines studied 
are not available, one may use an alternative 
expression for the intensity of a spectral line. 
This is, 

64r'y* = ea? 


cians 2. 6 
ny” (6) 


E a= ex — Xn k ] 
(4)? 


wo 3C* 





In this equation ¢ is the charge on the electron, 
a is the radius of the first Bohr orbit, 7 is the 
larger of the initial or final values of the quantum 
number in the two electron configurations from 
which the line arises, s is the “‘strength” of the 
spectral line, and p is the radial quantum in- 
tegral defined by 


ao 


yo J RuRnvr'dr, 
0 


where R,; and R,» are the radial wave func- 
tions. The factor p is a constant for all transi- 
tions between two electron configurations. If the 
spectral lines whose intensities are being com- 
pared belong to the same transition array, the 
factor p need not be included, since it is a con- 
stant for all the lines in this array. A large 
amount of data on theoretical s values is avail- 
able, and Petrie (see reference 4, p. 42) has 
shown how to obtain these factors from meas- 
ured laboratory line intensities. If the lines 
whose intensities are being compared belong to 
different transition arrays, then the p-factors 
must be included in the calculations. Menzel 
and Aller? have shown how these p-factors may 
be obtained in a number of cases. It is important 
to remember that the preceding equations apply 
only to an assemblage of particles in thermal 
equilibrium. Hence, if we use these equations 
with the intensities of lines arising from different 


5 L. Goldberg, Astrophys. J. 89, 673 (1939). 
6 W. Petrie, J. Roy. Ast. Soc. Can. 38, 137 (1944). 


7D. H. Menzel and L. H. Aller, Astrophys. J. 94, 436 
(1941). 
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TEMPERATURES DERIVED FROM EMISSION SPECTRA 


energy levels, we obtain information on the de- 
gree of departure from a condition of thermal 
equilibrium. 

The appearance in a spectrum of emission 
lines from the same element in two or more 
stages of ionization enables one to determine 
what is called the “ionization temperature.’”’ The 
number of neutral atoms of an element in a given 
energy level r may be obtained from one form of 
the Boltzmann expression, 


Nor Wor 


= —Xo0r kT e 
M Bot) exp(— xor/kT) 


(7) 
In this equation No, is the number of neutral 
atoms in level 7, No is the total number of neu- 
tral atoms. Divide (3) by (7) and we have 


Mey 2By(T) (24m)*!?(kT) 5? 


Nor Wor he 


Xexp(—(x—xor)/kT). (8) 


The Boltzmann equation may also be applied to 
the ionized atoms. The number of ionized atoms 
N,, in level s is 


Nie 
M 


is 


= exp(—x10/RT). 


B,(T) ” 


Forming the product of Eqs. (8) and (9), we have 
wis (2am)3/2(RT) 5/2 
P= SEE 
hs 
Xexp((—x+x0r—x1s)/RT). 


The relative intensities of the lines arising from 
these levels is 


N1,A 1s/V 1. 


NorA orlVor 


Aihvis 2wig (2am)?/2(RT) 5/2 
A orhvor Pwor hé 


seul ~ 2400 — va AT. 


The electron pressure in dynes cm~ is nkt, n 
being the number of free electrons per cm’. This 
method has been used to compute temperatures 
in stellar atmospheres. The ionization tempera- 
ture then, refers to the absorbed radiation or 
colliding particles which produce the observed 
equilibrium between ionized and neutral atoms. 
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Another quantity, the “electron temperature,” 
may be obtained from the continuous emission 
that in certain cases appears at the limit of a 
spectral series. Consider the origin of this con- 
tinuous emission. At high temperatures ioniza- 
tion and recapture processes are active. That is, 
atoms may absorb sufficient energy to lose one 
or more electrons. However, the resulting ions 
recombine with free electrons and become neu- 
tral atoms once again. For a given temperature, 
the Saha equation expresses the relative num- 
bers of neutral and ionized atoms. Since the 
energies of free electrons are unquantized, the 
radiation resulting from the recapture process 
will form a continuous spectrum. If we consider 
the continuous emission at the limit of the hy- 
drogen Balmer series, this emission being pro- 
duced by the capture of an electron on the level 


n=2, the frequency of the emitted quantum is 
given by 


hy =3mv?+thRe, 


where 3mv’ is the kinetic energy of the free elec- 
tron before capture, and (4)hRc the energy of 
level 2 measured from the ionization level, R 
being the Rydberg constant and c the velocity 
of light. The intensity distribution in this con- 
tinuous spectrum depends upon the velocity 
distribution of the free electrons; hence measure- 
ments of the spectral gradient in this continuum 
determine the electron temperature, that is, the 
temperature which produces the observed ve- 
locity distribution among the free electrons. 
Menzel and Cillié® have used this method to 
calculate the temperature of the solar chromo- 
sphere. 

In general, effective, color, excitation, ioniza- 
tion, and electron temperatures will differ for an 
assemblage of particles. Each temperature de- 
scribes a different phenomenon, and each con- 
tributes to our information on the physical 
conditions of the assemblage. Electron tempera- 
tures in certain cases may be less than ioniza- 
tion temperatures. The latter are determined 
from the relative intensities of lines from the 
ionized and neutral atom. After the ionization 
takes place the free electron may experience 
inelastic collisions with other particles and give 


8D. H. Menzel and G, G, Cillié, Astrophys. J, 85, 88 
(1937), 





382 Ww. J. 


up part of its energy. Hence, the energy emitted 
during the recapture process will be characteristic 
of a lowered temperature. If the upper levels of 
an ion in the nth stage of ionization are popu- 
lated by the recapture process, then the excita- 
tion temperature calculated from the spectrum 
of this ion will be much the same as the elec- 
tron temperature of the free electrons captured 
by the (m+1)st ionization stage of the ion. In 
other words the excitation temperature from the 
spectrum of O II will be close to the electron 
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temperature of the free electrons captured by 
the O III ion. 

We have pointed out that for an assemblage 
of atoms, molecules, and ions in thermal equi- 
librium all methods of computing temperature 
give the same result. However, thermal equi- 
librium conditions are not found in the earth’s 
atmosphere and in stellar atmospheres, and the 
choice of the various spectroscopic temperatures 
depends upon the physical phenomenon which 
one wishes to describe. 


An Experiment with a Nonlinear Circuit 


W. J. CUNNINGHAM 
Department of Electrical Engineering, Yale University, New Haven, Connecticut 


ONSIDERABLE time is commonly devoted 
to laboratory work with linear electrical 
circuits in courses of instruction dealing with 
circuit theory. Nonlinear circuits are often dis- 
missed with a few comments concerning the diffi- 
culty of analyzing them, even though many 
circuits of interest actually become nonlinear 
when large signals are applied. A relatively 
simple experiment involving forced oscillations 
in a nonlinear resonant circuit is described here. 
The experiment can be performed in a reasonable 
time with equipment usually available in a 
circuits laboratory. The nonlinear element of 
the circuit is an inductor with an iron core in 
which the magnetic flux is not a linear function 
of the current. The phenomena involved and the 
method of analysis are applicable to other 
resonant systems containing a nonlinear element, 
such as mechanical vibrating systems! or acous- 
tical devices.2 The method of analysis is stand- 
ard and is described elsewhere.* 

The circuit under study is shown in Fig. 1, 
and consists of an oscillator supplying a sinu- 
soidal voltage of variable frequency to the series 
combination of the iron-core inductor and a 


1S. Timoshenko, Vibration problems in engineering (Van 
Nostrand, 1937), p. 137. 
2H 


. F. Olson, Elements of acoustical engineering (Van 
Nostrand, 1947), p. 168. 


3L. A. Pipes, Applied mathematics for engineers and 
physicists (McGraw-Hill, 1946), p. 600, 
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fixed capacitor. The inductor must be chosen 
carefully, since the success of the experiment 
depends upon it. The core of the inductor should 
be of such design that it can be saturated mag- 
netically with a relatively small current in the 
coil, so as to avoid the necessity of large power 
input to the circuit. At the same time, the self- 
inductance of the coil should not exceed a few 
henries, and the ohmic resistance of the winding 
should be small. Data given here were obtained 
using the primary winding of a Western Electric 
transformer, type D-163413, which is designed 
to operate at 400 cycle sec with an input volt- 
age of 115 volts. The iron core is relatively small 
in cross section. Ordinary 60 cycle sec—! power 
transformers are not suitable unless a con- 
siderable amount of the core material is re- 
moved. Several small audio inductors have been 
found to function successfully in the circuit. The 
value of the capacitor should be such that the 
resonant frequency of the circuit is not too high, 
and the magnetic flux is therefore large. Volt- 
meters of high impedance are connected across 
the oscillator and across the inductor. An oscillo- 
scope connected across the inductor is also use- 
ful, since it allows observation of the wave form 
of the voltage at this point. 

If the driving voltage is maintained constant 
in amplitude, and its frequency is raised steadily 
from a low value, the rms voltage across the 
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EXPERIMENT WITH A NONLINEAR CIRCUIT 


inductor varies as shown by the solid curve of 
Fig. 2. Circuit parameters used in obtaining this 
figure are Eoy=15 volts rms, C=1yf, and the 
inductor described previously. The voltage Ex, 
rises with the frequency until, at a certain critical 
frequency (about 230 cycle sec for the figure), 
the voltage falls abruptly to a much smaller 
value, approaching the driving voltage asymp- 
totically. If the frequency is now lowered 
steadily, again maintaining constant amplitude, 
the voltage across the inductor increases only 
slightly until another critical frequency (about 
135 cycle sec!) is reached. At this frequency, 
the voltage jumps suddenly to the larger value 
of the curve obtained with increasing frequency. 
Between the two critical frequencies (indicated 
by the black squares of the figure) there is a 
range in which the voltage across the inductor 
may have either of two values, depending upon 
the past history of the circuit. Theoretically, a 
third value also exists in this range, but it is 
unstable and cannot be observed physically. 
The wave form of the voltage across the inductor 
becomes square-topped along the upper branch 
of the curve, for large values of this voltage. 
This distortion from sinusoidal wave form indi- 
cates the presence of a third harmonic com- 
ponent. 

A qualitative explanation of the behavior just 
described may be based on the observation that 
the effective inductance of the coil becomes 
smaller when the amplitude of the current in it 
is increased. As the frequency of the driving 
voltage is raised from a low value, resonance is 
approached. The resulting drop in the impedance 
of the circuit allows an increase in current, giving 
a corresponding decrease in inductance and a 
rise in the resonant frequency. This process is 
cumulative, with resonance continually being 
approached but never quite reached. Eventually 
losses in the system become so large that the 
high current cannot be maintained, and the drop 
in amplitude occurs. The new current is based 
on the original large inductance. If the frequency 


Fic. 1. Resonant circuit used in experiment. 
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Fic. 2. Rms voltage across inductor of experimental circuit 


for constant driving voltage, Eo=15 volts. 


is now reduced, resonance is approached once 
more, giving the rise in current and drop in 
inductance. Soon an unstable point is reached, 
such that a further rise in current would cause 
the inductance to become small enough to give 
a resonant frequency higher than the driving 
frequency, and the current jumps to the large 
value attained previously. The behavior of the 
circuit may be predicted by the approximate 
analysis which follows. 

It will be assumed that the resistance of the 
circuit is small and can be neglected. This as- 
sumption simplifies the argument considerably, 
since the actual resistance of the circuit would 
be a nonlinear quantity itself. Neglect of the 
resistance prevents the location of the frequency 
of transition from the upper branch to the lower 
branch, however. Approximately, the ‘relation 
between the current 7 in the inductor, and the 
magnetic flux linkage ¢ is 

4=$/Ly+a¢'. (1) 
Hysteresis effects are neglected in writing this 
equation. The self-inductance would be Lo if 
the coil were linear, and a is a coefficient deter- 
mined by the amount of nonlinearity. 

The voltage equation for the circuit of Fig. 1 
may be written as 


do /at+1/C f idt=E, sinwt. (2) 
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Fic. 3. Circuit for determining magnetization curve 
for inductor. 


The successive terms represent the voltages 
across the inductor and the capacitor, and the 
driving voltage. If Eq. (1) is substituted into 
the derivative with respect to time of Eq. (2), 
the result is 


d*¢/d?+¢4/LoC+a¢?/C=Ew coswt. (3) 


The solution of this equation cannot be found 
readily. However, experimental observations in- 
dicate that the voltage across the coil (and thus, 
the flux in it) varies with time in approximately 
sinusoidal fashion. Furthermore, if resistance is 
neglected, the driving voltage and that across 
the coil must either be in phase or in phase 
opposition. Therefore, a first approximation to 
the steady-state solution might be 


$= ¢1 coset. (4) 
Substitution of this solution into Eq. (3), to- 


04 
? 


weber turns 
03 


gether with the relations, wo?=1/LoC, b=a/C, 


and cos*wt = 3 coswt++4 cos3wt, gives 


(—wo1 + w0?1 + 2661) coswt +4541? cos3wt 
=Ew coswt. (5) 


If the inductor is nonlinear so that } is not 
zero, no value of ¢; will satisfy this equation. 
The difficulty arises from the fact that Eq. (4) 
is not a correct solution for Eq. (3). However, it 
is reasonable to require that the fundamental 
oscillation satisfy Eq. (5) and to disregard the 
third harmonic term for the present. Such a 
procedure yields the relation 


3613 = (w? — wo?) 1+ Ew. (6) 


This cubic equation relates the amplitude, ¢y, 
of the flux wave to the circuit parameters and 
the driving voltage. Since the instantaneous 
voltage across the inductor is e,=d¢/dt, the 
rms voltage at this point is 


E,=2-%e|¢1|. (7) 


In order to make use of Eq. (6) for predicting 
the performance of the circuit, the constants 
for Eq. (1) are needed. These constants may be 
found from a magnetization curve for the in- 
ductor. Such a curve may be determined experi- 


Fic. 4. Magnetization curve 
for inductor of experimental 
circuit. 
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mentally by the circuit of Fig. 3. The inductor 
is connected in series with a resistor, Re, of 
relatively small value, and an oscillator of ad- 
justable output voltage. The frequency is un- 
important but should be fairly low. An inte- 
grating circuit, Ri and Ci, is connected across 
the inductor. Typical values are Ri=0.25 meg, 
Ci:=0.5uf, Re=1000 ohms, and f=45 cycle 
sec". The horizontal and vertical amplifiers of 
an oscilloscope are connected as shown in the 
figure. The horizontal deflection on the oscillo- 
scope is proportional to the current, 7, and the 


vertical deflection is proportional to the flux 


linkage, ¢. It is necessary that the oscilloscope 
be calibrated so that the peak voltage required 
at the input terminals to produce a given peak 
deflection is known. Then the following relations 
apply: 

4=en/Ro, (8) 


6= f extt=RiCeer, (9) 


where ey and ey are instantaneous voltages at 
the input terminals of the horizontal and vertical 
amplifiers, respectively. 

The figure on the screen of the oscilloscope 
will be a hysteresis loop for the inductor. The 
loop will change in size and shape as the driving 
voltage from the oscillator, and the current in 
the coil, is changed. The tips of the successive 
loops will trace the desired magnetization curve. 
The curves for the D-163413 inductor are shown 
in Fig. 4. The complete magnetization curve is 
given by the solid line, while portions of several 
hysteresis loops are shown in the lower half of 
the figure. An approximation to the magnetiza- 
tion curve, as obtained from Eq. (1), is plotted 
in the first quadrant. The equation is too simple 
to allow the initial sharp curvature near the 
origin to be fitted, but a fairly good fit with the 
upper portion of the curve is possible. It is best 
to choose constant Lo so the approximate curve 
has a slope at the origin about the same as the 
maximum slope of the actual curve. Constant a 
is then chosen so as to fit the knee of the actual 
curve fairly well. The equation used in Fig. 4 is 


i=$/7.2+125¢%, (10) 


where 7 is measured in amperes and ¢ in weber 
turns. 
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Fic. 5. Graphical solution for amplitude of flux linkage 
for experimental circuit. 


A substitution into Eq. (6) of the known quan- 
tities for the experimental circuit gives 


9.4 X 107613 = (w? — 1.39 X10°)¢14+21.20. (11) 


The flux and voltage in this equation are peak 
values. A graphical solution of the equation is 
shown in Fig. 5. The solid curve of the figure is 
the cubical parabola representing the left side of 
the equation. The dotted straight lines represent 
the right side of the equation for various fre- 
quencies, as indicated. The intersection of these 
lines with the vertical axis is determined by the 
second term of the right side, while the slope of 
the line is determined by the first term. Values 
of ¢1, found at the intersections of the solid 
curve and the dotted lines, are solutions of Eq. 
(11). For relatively low frequencies, only a 
single intersection occurs. Above a certain fre- 
quency, three intersections occur. The dotted 
line for 140 cycle sec~! shows two intersections, 
and a third would occur if the figure were ex- 
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Fic. 6. Curves analogous to those of Figs. 2, 4, and 5, but using a more accurate approximation 
to the measured magnetization curve. 


tended sufficiently. The rms voltage across the 
inductor at any frequency is given by Eq. (7), 
where the appropriate value of ¢; is taken from 
Fig. 5. Points corresponding to values of Ex, de- 
termined in this way are plotted as circles in 
Fig. 2. The calculated curve drawn through these 
circles is seen to agree fairly well with the 
measured curve. The portion of the right-hand 
branch of the curve having a positive slope is 
unstable. Therefore, the jump in amplitude at 
135 cycle sec! as the frequency is lowered is 
predicted. Since resistance has been neglected, 
the drop in amplitude as the frequency is raised 
cannot be predicted. 

The dotted curve of Fig. 2 is the calculated 
behavior if the inductance had the constant 
value, Lo. This curve may be determined from 
Fig. 5 by allowing the cubical parabola to de- 
generate to the horizontal axis. The curve for 
the linear circuit is very different from that for 
the nonlinear circuit. 

The calculated and measured voltages for the 
nonlinear circuit may be brought into better 
agreement at the lower frequencies by a further 
correction. Equation (3) and Eq. (4) may be 


combined to give 


d?¢/di? = (Ew - word os 3bd;°) coswt 
—+3b¢1' cos3wt. (12) 


A substitution of Eq. (6) gives 

d*¢/dé = —w’¢1 coswt —4b¢;' cos3wt. (13) 
Integration once with respect to time yields 
e,=d¢/dt = —w¢; sinwt — (bd13/12w) sin3wt. (14) 


The first term is the first approximation of the 
fundamental frequency, already used. The sec- 
ond term is a correction term having the third 
harmonic frequency. Its rms value is 


E13) = 27*bg13/12w. (15) 


The phase relations indicated by Eq. (14) are 
such as to give a square-topped wave. Approxi- 
mately, the peak of ez is less than the peak of 
the fundamental component by the peak value 
of the third harmonic. The vacuum tube volt- 
meter used in measuring Ez, as plotted in Fig. 2, 
measured peak value although its calibration 
was in rms value. Therefore, its reading was too 
low if only the fundamental component is con- 
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sidered. Values of Ex:3), as found from Eq. (15), 
may be subtracted from the values of Ez, found 
from Fig. 5 and Eq. (7), to give more nearly the 
values indicated by the voltmeter. These cor- 
rected points are plotted as crosses in Fig. 2. 

The region of largest difference between 
measured and calculated curves of Fig. 2 is 
along the upper branch where the point of dis- 
continuity is approached. In this region the 
value of ¢: is so large (see Fig. 5 for 100 cycle 
sec!) that the approximate curve of Fig. 4 is 
no longer close to the actual magnetization curve. 
A better agreement between these curves may 
be obtained by the addition of a third term to 
Eq. (1), so that it becomes 


t= $/Lo+m¢*?+ng¢'. (16) 


The three coefficients of Eq. (16) allow it to fit 
the measured magnetization curve more closely 
than was possible with Eq. (1), but considerably 
more mathematical work is now required to find 
these coefficients. They must be chosen care- 
fully or a poor fit will result. Successive terms of 
the equation are large and of opposite algebraic 
sign, while their sum remains small. A poor 
choice of the coefficients causes the curve of 
Eq. (16) to cross the magnetization curve at 
certain chosen points, but to oscillate violently 
either side of the curve between those points. 
A fairly good fit between the curves is shown in 
Fig. 6(a), where the equation is 


4=/4—10°o*+ 10°9°. (17) 


It is not difficult to show that, following the 
previous analysis, the amplitude of the flux 
linkage is now given by 


3m¢1*/4C+5ng15/8C = (w*— wo?) di tEw. (18) 


The solution for ¢1 may be carried out graphi- 
cally as before. This process is shown in Fig. 6(c), 
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and the resulting curve of voltage across the 
inductor is shown in Fig. 6(0). No additional 
correction has been applied, but the calculated 
and measured curves here are seen to agree 
more closely than do those of Fig. 2. 

The behavior of the circuit of Fig. 1 is actually 
somewhat more complicated if small driving 
voltages are considered. With very small driving 
voltage, the inductor is nearly linear and a 
simple resonance peak is observed. However, 
the effective value of the inductance at very 
small currents is considerably less than the con- 
stant Lo of Eq. (10) or Eq. (17). In fact, this 
inductance is given by the slope of the magnetiza- 
tion curve of Fig. 4 at the origin, and is less 
than one henry for the coil used here. The 
resonant frequency with small currents is there- 
fore much higher than the wo of the analysis. 
As the driving voltage is increased, operation 
swings up the rising portion of the magnetiza- 
tion curve, and the effective inductance becomes 
larger, thereby lowering the resonant frequency. 
Still further increase in the driving voltage 
causes operation over the knee of the magnetiza- 
tion curve, and the phenomena previously de- 
scribed take place. An analysis of the complete 
behavior would require still more terms in Eq. 
(1), and the entire procedure would be much 
more complicated. 

In spite of the simplifying assumptions used 
in the experiment, it serves to illustrate the 
measurement and analysis of a typical non- 
linear resonant system and leads to results which 
are fairly good. The techniques of measuring 
the magnetization curve for the iron-cored coil, 
and the graphical solution of the equation for 
the amplitude, are interesting tools for other 
purposes. The performance of the experiment 
and the necessary calculations can be carried 
out in a reasonable length of time. 


The instinctive weapon of the bourgeois, the reading and writing class, ts the book. Whenever 
he is in a state of danger or hope he starts furiously to learn.—WILL1AM BOLITHO. 



















































































L. E. GRINTER 
Illinois Institute of Technology, Chicago, Illinois 


T seems to me that anyone with the temerity 
to suggest improvements in the handling of 
graduate study had best start by examining his 
own qualifications. This is especially true because 
the research professor and the dean are likely to 
see the over-all problems of graduate study from 
very different perspectives. Administration repre- 
senting organization, and teaching, which sym- 
bolizes inspiration, are more widely separated 
in graduate than in undergraduate education. 
If these natural divergences have meshed with- 
out conflict in my own thinking it is, perhaps, 
because I have had the opportunity to complete 
the cycle of graduate teaching, graduate ad- 
ministration and graduate teaching again. It 
has been no small contribution to my educa- 
tion to return during the past two years to 
teaching and research where it has been possible 
in an impersonal way to observe the workings of 
one of the larger graduate schools of engineering. 
It is unfortunate that so few presidents, deans 
or professors have the opportunity to complete 
this cycle of observation. It is a sobering and 
calming experience. 

If a re-inspection of graduate study is to prove 
useful, it should search out the weaknesses of 
graduate education and strengthen them. But 
weakness or strength depends upon the end being 
sought. A practical viewpoint can be a weakness 
if the objective of graduate study is mathe- 
matical facility, but it may be a mark of strength 
in a curriculum designed to forward a construc- 
tive art. And surely it must be agreed that the 
broad sweep of the physical sciences and their 
applications through engineering include the en- 
tire range from abstract theory to the most ap- 
plied type of design. Of course, a part of these 
studies belongs to undergraduate education and 
another part may be found only in design offices. 
In any case, graduate study cannot confine itself 
wholly to the mathematical or theoretical studies 
even including experiment. It must also make 
provision for the less routine or higher aspects 
of design investigations. 

When the objectives of graduate study in any 


First-Line Problems of Graduate Study 
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particular department are being settled, the 
greatest influence thereon should be exerted by 
the interests of the teachers. There is no reason 
why a department of physics, for example, must 
teach a graduate course in spectroscopy. The 
decision should depend upon whether there is a 
teacher on the staff who has a special interest and 
ability to teach such a course; and my criterion 
would be whether research in spectroscopy was 
actually under way in that department. I find 
no difficulty in commending a graduate school 
in civil engineering devoted exclusively to hy- 
draulics without the offering of a single course in 
structures, or vice versa. Students should select 
a graduate school on the basis of its reputation 
for highly specialized studies. 

You will see then that I am not in sympathy 
with the idea that graduate study may be just a 
fifth year of education in a broad field. There is 
much in common in my thinking about a fifth 
year and a fifth wheel. We start a group of 
freshmen that includes within it (a) future scien- 
tists, or specialists, (b) future technological sales- 
men, promotors and managers, (c) future con- 
struction and production men, and (d) many 
young men who think they would rather study 
science or technology than liberal arts, commerce 
or agriculture. For three of these four groups, 
four years of study is the equivalent of four 
wheels on a cart—it serves all possible purposes. 
For some, three years or three wheels provide an 
adequate educational vehicle. For others a two 
wheeled education serves just as a two wheeled 
cart will transport necessities. Even the man 
with one year of college education has at least 
the aid of an educational wheelbarrow which is 
far superior to loading one’s back! 

No doubt some one will point out that the 
fifth wheel was an important invention when 
added to the wagon since it permitted the opera- 
tor to guide the vehicle. And this too seems to 
me an apt analogy in education. It is worse than 
wasted to add a fifth wheel just anywhere on the 
vehicle. But when scientifically placed in a hori- 
zontal plane it adds a new quality of direction or 
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guidance, the perfect symbol of vision or leader- 
ship in science and technology. Leadership in 
science develops only from deep knowledge 
gradually broadened through experience. A gen- 
eral fifth year of broadening education uni- 
versally required in all technological curricula 
(and commonly preceeding the junior year) 
might be likened to putting a five wheel cart in 
front of a horse without a driver. Nevertheless, 
if adopted only here and there, such a fifth year 
might meet the needs of limited numbers of 
students who would then select such curricula 
because of their special interests. I criticize only 
the universality of the idea as included in the 
common question ‘When will we come to five- 
year curricula in engineering and science?” 

In June 1945 the Journal of Engineering Educa- 
tion printed the Manual of Graduate Study in 
Engineering, a 37-page report of a three-year 
continuous study by a committee of experienced 
graduate educators. As I reread this manual it 
seems to me that it is as applicable to any field 
of science as it is to engineering because any 
committee of competent educators must in- 
evitably come to about the same conclusions. 
Without being didactic it tries to outline good 
procedures and desirable objectives concerning 
the matters of admission requirements, degree 
requirements, major and minor studies, the 
thesis, language requirements, mathematics, non- 
technical studies, undergraduate courses for 
graduate credit, required examinations, and so 
forth. In review I still find that it performs a 
very useful service. 

Nevertheless, I also find that the manual over- 
emphasizes the techniques or the detailed or- 
ganization of graduate study in relation to the 
less tangible but more significant factors. This is 
natural because it is easy to be definite and to 
write definitely about the need for written and 
oral examinations while it is difficult to think 
about and still more difficult to write about the 
intangible quality of inspiration of students by 
the faculty. Nevertheless, the latter is important, 
and the kind of examination used is trivial by 
comparison. If the examination becomes in- 
effective, as I have observed occurring in my 
own institution, it can be made effective by 
changing its form, by transferring its position in 
the curriculum, or by using a different system 
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of selecting the examination committee. But 
what is to be done if students are not being in- 
spired by the faculty in a whole department or 
school? How are new faculty members to be 
selected? Is mediocrity to be permitted to per- 
petuate mediocrity as it will, inevitably? Are 
young staff members to be encouraged, stimu- 
lated and rewarded for accomplishments in re- 
search? Are the few, the very few, great graduate 
teachers to be rewarded in proportion to their 
distinction or are our institutions to gamble on 
their loss in order to maintain a balance in 
institutional salaries and teaching loads? In the 
answers to such questions lie the more significant 
influences upon graduate study. Such answers 
lead to one institution becoming great in the 
national sense while another remains effective 
merely in its own locality. 

If an institution is interested exclusively in 
undergraduate education, or nearly so, it can 
avoid solving many of these embarrassing prob- 
lems. But the graduate institution must face 
these issues squarely. Administrative officers 
carry a large share of the responsibility. A great 
president will attract great scientists to his in- 
stitution and then he will contrive to have them 
shine as brilliantly as himself. An effective dean 
will be known by the reputations of his key 
teachers. Beware of any institution which is 
known to have a brilliant president or dean but 
an unknown faculty. There are even department 
heads who avoid the competition of brilliant 
associates; but the best administrators seek such 
competition as the life blood of an institution. 

After a distinguished graduate teacher and 
researcher or a promising younger scientist is 
attracted to the institution, what further effort 
should be made to assure his continued research 
productivity, without which his teaching will 
soon become stale and eventually will stop being 
fruitful? Here is the place where education has 
the least reason to be congratulated. I recall 
reading that the investment in tools for a factory 
worker averages about $6000 and _ reaches 
$50,000 in some industries. If such an invest- 
ment is needed to make an average workman 
productive, is it not reasonable to expect a 
scientist to require an investment many times 
as great? The institution that hires a graduate 
professor and assumes that he can be made pro- 
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ductive by providing him with a private office 
and little else is due for disappointment. 

The theoretical analyst usually needs one or 
more assistants who are more than inexperienced 
students. He may need calculating machines of 
a simple or complex nature. He always needs 
secretarial service that is not merely the right to 
borrow the help of another person’s secretary. 
With less than these his productivity will be 
gradually strangled. Yet how many universities 
really provide this minimum standard? The 
experimentalist requires a laboratory and a shop 
with laboratory assistants and a mechanician. 
Everyone has seen professors, young and old, 
making their own equipment at the expense of 
time that cannot be purchased. To the industrial 
research director this self-service merely looks so 
foolish that he doubts its existence. Industrial 
research laboratories are overstaffed with hands 
and understaffed with brilliant minds. Relatively 
speaking, the converse is true in educational 
research. 

How many educational institutions provide a 
service whereby a professor can get high quality 
drafting done without taking time to train the 
draftsman himself? If a professor does not have 
full time need for a mechanic, an electrician or 
an instrument maker, can he obtain such ser- 


LLOYD WILLIAM TAYLOR 


vices from a central service bureau or must he 
beg the loan of a workman from an associate 
who runs a laboratory, or from the buildings and 
grounds department? Industrial people cannot 
understand why we show such lack of foresight 
in improving our research output. 

The war years produced at least one good re- 
sult in educational research. The availability of 
government funds was used for the first time in 
many institutions to provide collateral services 
to research workers. Postwar contracts have 
helped to keep such services in existence. Also, 
the exchange of workers between education, in- 
dustry and government laboratories served to 
educate college administrators and led research 
workers to insist upon more aid in research. 
Nevertheless, this lack is still predominant, and 
it reflects little credit upon the efficiency of 
educational operations. 

It is such broad questions that lie at the heart 
of the problem of producing effective leadership 
in graduate study and research. We are always 
short of brilliant minds. Let us try to make those 
that education is lucky enough to have as ef- 
fective as possible by providing them with the 


necessary tools and assistance. Only in this way 
can the hours be freed that are so badly needed 
for real direction of graduate study and research. 


Lloyd William Taylor 


We regret to have to record the tragic death of Lloyd 
William Taylor, Professor of Physics at Oberlin College, 
who was killed in a mountaineering accident on Mt. 
St. Helens, Washington, on Sunday, August 8, 1948. 
To commemorate in appropriate fashion his notable 
contributions to physics and to the American Associa- 
tion of Physics Teachers, a Lloyd William Taylor 
Memorial number of the American Journal of Physics 
is being planned for the early part of 1949.—T. H. O. 





A Null-Deflection Magnetometer with Electromagnetic Control 


A. R. INGLES 
School of Physics, College of Technology, Leicester, England 


HE couple exerted on a coil with its plane 

in a magnetic field is directly proportional 
to the current in it. If an external torque be 
applied to such a coil, causing it to turn, it 
follows that the current required in the coil to 
restore it to its original undeflected position is 
proportional to this external torque. Hence it 
should be possible to convert a center-zero 
moving-coil detector galvanometer, such as is 
usually employed in a Wheatstone bridge cir- 
cuit, so that it could be used either for measuring 
small forces or torques or as a null-deflection 
magnetometer. 

The purpose of this article is to explain how 
the latter aim has been accomplished and to 
describe some experiments carried out with the 
magnetometer. It is hoped that it may serve to 
stimulate the interest of those teachers and 
students who have the desire and the time to 
make and use apparatus of a somewhat uncon- 
ventional kind, for in so doing much valuable 
knowledge and technical skill can be acquired. 
The particular instrument described here is 
somewhat fragile and rather large, but there is 
no reason why improvement could not be made 
in these respects by the choice of other and per- 
haps more suitable components and materials. 
The magnetometer has a particular advantage 
in that its operation is independent of the direc- 
tion of external magnetic fields, provided they 
do not vary during the time of the experiment. 


Details of the Magnetometer 


The moving-coil instrument selected for use 
had a pointer with a center zero scale. It was 6.5 
cm high and its 10-cm diameter dial was cali- 
brated 5-0-5 ma. It was an old instrument and 
its maker could not be identified. 

It was first tested to see whether the currents 
needed to restore the pointer to its zero mark 
were directly proportional to the loads applied 
to the end of a thin aluminum rod fixed by 
paraffin wax to the coil mechanism so that it 
projected beyond the case of the instrument, 
opposite to and in line with the pointer. The 


weight of the aluminum rod was counterpoised 
by fixing lead shot to the pointer by paraffin 
wax. The heated blade of a small screwdriver 
was found to be a convenient tool for manipulat- 
ing the wax. The instrument was fixed to a board 
so that its pointer moved in a vertical plane. 

For calibration purposes, six riders were cut 
from copper wire and weighed. Each was applied 
as a load at a fixed point on the aluminum rod 
and the current required to bring the pointer 
back to the no-load mark was recorded. The 
values of load and current given in Table I and 
shown graphically in Fig. 1, show clearly that 
these quantities are directly proportional. 

In this form the instrument could well be 
used to measure forces ranging from 1 to 23 
dynes. There is no reason why a more sensitive 
balance could not be made by choosing a more 
suitable coil movement. 

The conversion of the chosen movement into 
the magnetometer used in these experiments was 
next carried out. The arm / replaced the thin 
aluminum rod already mentioned. It was made 
by drawing out a glass tube and bending it to 
the shape shown in Fig. 2. It was fastened to the 
movement C by means of paraffin wax as before. 

Two small Alcomax 2 magnets! M,, approxi- 
mately 5X11 mm, were mounted at the end 
of arm h, also by wax. A concave mirror m of 
radius 100 cm was fixed to the arm over the 
axis of rotation of the movement. 

The load added by these parts was counter- 
poised by means of lead shot L stuck onto the 
pointer P by wax. The final adjustment for 


TABLE I. Restoring current and mass of rider. 


Rider Restoring 


ass current 
(mg) (ma) 


4.0 0.84 
1 1.61 
11.5 2.43 
15.1 3.20 
19.0 4.03 
22.7 4.81 


1 Supplied by Messrs. W. Jessop & Sons of Sheffield. 
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balance was made by moving the aluminum rod 
A in or out of the end of the hollow glass arm h, 
to which it was finally fixed also by wax. The 
balance thus obtained was correct for both ver- 
tical and horizontal movements of the moving 
system. 

The arm h/ was about 32 cm long and was 
made as light as possible yet rigid enough to 
prevent sag under the load of the magnets M, 
and also to prevent whip when subjected to the 
normal vibrations of a laboratory bench. 

Figure 3 shows in plan and elevation how the 
magnetometer was mounted in a small cabinet 
with a hinged lid, fitted with a glass front and 
side. The cabinet was 23 cm wide, 27 cm high 
and 55 cm long. The box was essential to elimi- 
nate the effect of air currents. It should be noted 
that the movement was now in a horizontal 
plane. 





Fic. 2. Method of attaching the small magnets to 
the coil movement. 


Fic. 1. Linear relationship between the current required to restore the moving coil to its undeflected position 
and the load applied to a fixed point on the arm attached to the coil. 






The box was fitted with plane mirrors X and 
Y to enable the deflecting magnets to be lined 
up correctly on the axis of the magnets M). A 
hole H was cut in the lid of the box so that be- 
fore each experiment the no-deflection position 
of M, could be checked and reset if necessary by 
the zero control on the movement; however, no 
such adjustment was found to be necessary. The 
distance between the center of M, and the glass 
front was accurately measured. 

The restoring current required to bring the 
deflected magnets M, back to the no-deflection 
mark was controlled and measured with the 
circuit given in Fig. 4, the terminals of which are 
connected to those of Fig. 3. In Fig. 4, Gis a 
Cambridge Versatile galvanometer of resistance 
500 ohm, reading to 120 mv, and shunted by a 
resistance box S; R is another resistance box and 
C is a commutator; E is a 2-v storage cell and 
R, a 19-ohm, 2.8-amp rheostat used as a poten- 
tial divider. The resistances of S and R, were 
adjusted so that for a given experiment the cur- 
rent was entirely controlled by the movement of 
the slider of Ro. 

The damping of the magnetometer was such 
as to cause no delay in taking readings and when 
the key K was closed it was considerable. 
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A NULL-DEFLECTION 


Experiments Performed 


(1) Force between a very small magnet and a 
short magnet placed on a common axis 


Theory of experiment.—lf the magnetic field 
strength at P (Fig. 5) be H, then the force on the 
very small magnet of length dx, pole strength 
m,; and moment M, is given by 


(dH/dx)dx; "m= (dH/dx)M,, 
and thus 


(dH/dx)Myh=bi, (1) 


where 6 is the couple per unit current acting on 
the coil of the movement which is carrying cur- 
rent ~ to maintain no deflection. 


For a magnet of moment M and magnetic 
length 21, 


H=2Mx/(x?—P)? 


and 
dH/dx = —2M(3x?+P)/(*?—P)3. 


If the magnet is small, so that P is negligible 
compared with x?, 


dH/dx = —6M/x*, 


whence the force on M,; due to M is 6MM,/x*. 
Thus 
6M M,h/x*=0i, 
whence 
1x‘ = const. (2) 
and 
Int+4Inx =const. (3) 


Experimental details —The Alcomax 2 magnet 
M (Fig. 5) had dimensions 4X1X1 cm and was 
supported on a platform made by mounting a 
meter stick on wooden supports at such a height 
that it lay with its axis along the axis of the 
small Alcomax 2 magnets M, on the mag- 
netometer. The correct positions of the platform 
and the magnet M were found by means of the 
images reflected in the glass front and in the 
vertical mirror on the remote side of the mag- 
netometer case. 

Values of restoring current 7 were measured 
for values of x between 13 and 19 cm, and for 
both polarities of magnet M. The experimental 


values of ix* differ from the mean by 3 percent 
or less. 


MAGNETOMETER 


Fic. 3. Plan and elevation of arrangement used in 
mounting the magnetometer in its cabinet. 


(2) Force between a very small magnet and a 
single pole placed on the axis of the magnet 


Theory of Experiment.—The field strength at 
P (Fig. 6), due to the Robison magnet NS of 
pole strength m, is 


m™ m 


—— cos@ 
x? (x? op L?) 


m mx 


x (x2-L2)8/2" 


L?—2x? 
m\| ————— |. 
ae 
For values of x which are small compared with 
L the second term may in the right-hand member 


be neglected in comparison with the first, whence 
the force on M, due to the single pole m is 


Fic. 4. Circuit used for controlling and measuring 
the restoring current. 
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Fic. 5. Arrangement of the deflecting magnet and the 
magnetometer in Experiment (1). 


Fic. 6. Arrangement of the Robison magnet and the 
magnetometer in Experiment (2). 


2mM,/x*; thus 


2m 
—M,h=bi 
x3 


1x* = const. 
Int+3 Inx =const. (5) 

Experimental details—The Robison magnet 
NS (Fig. 6) was 44 cm long and had pole m 
supported on the wooden platform already men- 
tioned; to obtain a constant and large value of 
the strength of this pole the magnet was placed 
in a solenoid 41 cm long wound on glass and 
having two layers of turns of No. 24 S.W.G., 
double cotton covered copper wire wound 12 
turns to the centimeter. 

A current of 3.5 amp was found to be sufficient 
to saturate the steel of the magnet; it was re- 
duced to 0.5 amp during the experiment. It was 
also found that 0.5 amp produced no deflection 
on the magnetometer when the solenoid was 
used without the magnetic core, even with the 
smallest values of x. The values of x used lay 
between 10 and 15 cm. 


Experimental values of ix* were again found 
to be very nearly constant. 


(3) Magnetization curve and hysteresis cycle 


Experimental details —The specimen used was 
a knitting needle K (Fig. 7) to which was soldered 
a copper handle H. It was placed, after demag- 
netizing it in a separate alternating-current 
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circuit, in the solenoid LZ used in Exp. (2), with 
one end opposite the magnetometer magnets M). 

A 21-ohm, 4 amp rheostat was used as a 
potentiometer with an 18-v battery; the am- 
meter A had a range of 5 amp and the current 
was reversed by means of the commutator C. 

The effect of the solenoid was considered 
negligible since when it was carrying 4 amp it 
caused a deflection of only 2 divisions and no 
deflection at all with 1 amp as compared with a 
deflection of 120 divisions due to the needle K 
when it was saturated. 

Great care had to be taken in that part of the 
cycle when the polarity of K was opposite to 
that of M,, since the position of no deflection 
of M, was then one of unstable equilibrium. 

Readings of the solenoid current and the 
restoring current to bring the magnetometer 
back to its initial no-deflection position were 
taken over a complete cycle, the former being 
proportional to the magnetizing force and the 
latter to the intensity of magnetization. The re- 
sults are shown graphically in Fig. 8. 


(4) Variation in intensity of magnetization 
with temperature 


Experimental details —-A small bar magnet of 
dimensions 5.00.75 X0.6 cm was placed with 
its axis in line with that of the magnets M, of 
the magnetometer and was heated in the furnace 
F shown in Fig. 9. The furnace was made by re- 
placing the five burners in a standard tube fur- 
nace, 15 cm long, 10 cm wide and 18 cm high. 
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Fic. 7. Arrangement of the specimen, the magnetizing 
solenoid and the magnetometer together with the circuit 
used to control and measure the magnetizing current in 
Experiment (3). 
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A NULL-DEFLECTION MAGNETOMETER 


Fic. 8. Representation of the 
magnetization curve and hys- 
teresis cycle for the steel knit- 
ting needle used in Experiment 
(3). The restoring current is 
proportional to the intensity of 
magnetization and the magnetiz- 
ing current to the magnetizing 
force. 


by a single copper tube C, 34 cm long and 1.5 
cm in internal diameter, drilled at one end with 
four rows of six holes each to serve as jets after 
the manner of a gas poker; the tube was sup- 
ported in a clamp at the opposite end and gas 
was admitted by a bunsen burner B as shown 
and the gas flow was controlled by a screw 
clamp K. 

The end of C was plugged with asbestos wool 
A and the ends projecting through the firebrick 
of the furnace were passed through asbestos 
disks, thus keeping the heat at the end facing 


estoring Current (divs) 


the magnetometer from its glass front. The dis- 
tance x was about 17 cm. 

The temperature was controlled by the size 
of the flame jets and by putting on or taking 
off the three firebrick covers D; it was never 
necessary to make the flame jets very large. 

The intensity of magnetization of the steel of 
the magnet was proportional to the restoring 
current used with the magnetometer. The tem- 
perature was measured with an iron-constantan 
thermocouple XY connected in series with a 
resistance box R and a Cambridge Versatile 
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Fic. 9. Arrangement of the deflecting magnet, the gas 
furnace, the thermocouple and the magnetometer used in 
Experiment (4). 


galvanometer G of 10-ohm resistance and 2.4-mv 
range. 

Results are shown in Fig. 10 for two new 
magnets. Both magnets were made of the same 
steel but its exact nature was not known. 

Comments.—It is interesting to note the re- 
versal of polarity at W for the first magnet and 
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to attempt an explanation of it. (For the second 
magnet the temperature was reduced before 
the critical value was reached.) 

If it be assumed that for the steel of this mag- 
net the permeability increases rapidly to a 
maximum as the critical temperature is ap- 
proached and then falls quickly to unity at the 
critical temperature for small values of the mag- 
netizing force, as is the case for iron,? then the 
magnetism induced in the bar magnet by the 
magnets M, of the magnetometer might be 
sufficient to overcome that due to the earth’s 
field together with the permanent magnetism 
remaining after the heating to which the bar 
had already been subjected. 

If the temperature had been further increased, 
the curve beyond the point Y (Fig. 10) would 
presumably have followed the shape shown by 
the dotted line, ending by meeting the tempera- 
ture axis at Z, the ferromagnetic Curie point. 

It should be pointed out that the horizontal 
component of the earth’s field made an angle of 
about 65° with the axis of the bar magnet and 
was in such a direction as to cause it to retain 
its original polarity, the effective magnetizing 
force thus being Hy cos65°. 

On cooling from the point Y the reversed 
polarity acquired as suggested above would be 
retained but would decrease, assuming that the 
permeability decreased with fall in temperature 


Fic. 10. Variation of the intensity of magnetization with temperature for the two new steel magnets used in Experi- 
ment (4). The restoring current is proportional to the intensity of magnetization. +: Magnet 1, temperature increasing. 
©: Magnet 2, temperature decreasing. X : Magnet 3, temperature increasing. E]: Magnet 2, temperature decreasing. 


2R. M. Bozorth, Rev. Mod. Physics 19, 57 (1947). 
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as it does for iron in weak fields. Again if it be 
assumed that the field of the magnets M, in the 
region of the bar magnet was comparable with 
the effective field of the earth, Hy cos65°, then 
as these were in opposite directions the previous 
explanation of the reversal of polarity at W as 
due to induction is unsatisfactory. 

If, however, the bar magnet had been origi- 
nally magnetized with polarity opposite to that 
possessed by it at the beginning of the experi- 
ment, then on heating it might reverse its 
polarity when nearing the critical temperature 
owing to the persistence of some of the more 
stable domains, oriented in their original direc- 
tion, after the disappearance of the less stable 
domains oriented by the subsequent and re- 
versed magnetizing force. 

This explanation was suggested on reading of 
the problem of providing stability in the mag- 
netization of vessels after being subjected to 
the process of ‘‘wiping’’ or stroking with a mag- 
net or conductor carrying a current, developed 
in wartime for the purpose of changing the origi- 
nal magnetization in the steel hulls of ships as 
a means of protection against magnetic mines.* 

Although the former explanation is the more 
likely, as the magnet was a new one, it may be 
that the reversal of polarity in question was in 
some measure due to both causes. 


Other Possible Experiments 
(i) To compare the moments of two magnets.— 
From Eq. (1), 
(dH/dx) Myh=bi, 
where 
dH/dx = —-2M(3x?+P)/(x2—-P)*; 
3S. H. Ayliffe, J. Inst. Elec. Eng. 93, Part 1, 509 (1946). 
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era ocloei. a 


thus 
2M(3x*+P) h 


1- +1. 


(x* —F) b " 


Then for magnets of moments M: and M3 and 
lengths 2/, and 2/3, each at distance x from the 
magnetometer of moment M,, 


2M2(3x?+1,") (x? —13?)8 14 
(x# —1,?)8 


2M3(3x?-+13%) iz 
M2 1 (3x?+-13?) (x? —1,*)8 14 


Ms is (3x?-+12) (x—1g)? in 


for short magnets or magnets of equal length. 

(ii) To find the moment of a magnet.—It would 
first be necessary to find the constant of the 
magnetometer M,h/b shown in Eq. (6). 

If a flat coil of radius a and N turns carrying 
a current J amp were used, its axis being in 
line with the axis of the magnetometer magnets 
M, and its center at distance x from the center 
of the magnets M,, then the field due to this coil 
would be 


ete! 
~ 10(a2-+-x2)8/2" 
dH  6nxNI 


dx 10 


a’x 
(a2-+-x2) 5/2" 
and since, from Eq. (1), 

(dH/dx) M,(h/b) =7, 


the value of Mih/b could be determined, since 
both z and dH/dx are known. 

Once the value of this constant was established 
it could be used in Eq. (6) to find the value of M. 


the teaching of science to the young, so that 


when the whole population grows up there would be a far more general background of common 
sense, based on a knowledge of the real meaning of the scientific method of discovering truth. 
—Davip O. Woopsury, Beloved scientist, a life of EL1HU THOMSON. 





A Problem about Moving Charges 


Han Konc-Cui 
National Chekiang University, China 


T is generally agreed that electromagnetic 
phenomena can be best explained by Max- 
well’s theory of the electromagnetic field, in 
which electric and magnetic actions are con- 
sidered as propagated with finite velocity through 
a medium. In most problems, however, we are 
concerned with charges and magnets at rest or 
moving with steady small velocities. It happens 
that in such problems the assumption of an in- 
stantaneous “‘action at a distance” leads to the 
same results as an action propagated with finite 
velocity, and it is, moreover, a much simpler 
assumption to apply. As a result, many students 
fail, in their courses, to get a proper grasp of the 
role of the medium in Maxwell’s theory. 

The inadequacy of the theory of instantaneous 
actions can be most easily illustrated by prob- 
lems involving swiftly moving charges. One of 
the simplest of these problems is discussed in 
this note. 

Let two equal charges A and B, at a distance 
d apart, move with velocity v parallel to the x 
axis. From the point of view of an observer at 
rest, each charge will exert both an electrostatic 
and a magnetic force on the other. The magnetic 
force, computed on the basis of instantaneous 
action at a distance, will be 


F=H(ev/c) =ev?/d?c?. (1) 


Here, e is the charge in esu and H is the magnetic 
field in emu at one charge due to the other, as 
computed by the law of Biot and Savart. From 
the point of view of an observer moving with the 
charges, the magnetic force will be zero. It fol- 
lows that the resulting acceleration will be dif- 
ferent for different observers moving at constant 
relative velocity. This is, of course, a paradoxical 
conclusion. If v is restricted to values small com- 
pared with c, the magnetic force will be neg- 
ligible and the difficulty will not arise. 

The difficulty is avoided entirely if the forces 
and accelerations are computed on the basis of 
action propagated at finite velocity. To show 
this, let us find the components of the electric 
and magnetic fields at A and B, first from the 


point of view of an observer moving with the 
charges, and then, with the aid of a Lorentz 
transformation, from the point of view of an 
observer at rest. The former will be 


E'(A)z=E"(A).=0, E’'(A)y=—e/@; 
H'(A):=H"(A),=H'(A).=0; 

E'(B)z=E"(B).=0, E'(B)y=e/@; 
H'(B),=H'(B),=H'(B):=0. 


A Lorentz transformation applied to these yields 
for the components from the point of view of the 
observer at rest: 


E(A):=E(A).=0, E(A)y= ~(e/d’)(1— 62); 
H(A):=H(A)y=0, 
H(A).= — (Be/@)(1— 8-4; 
E(B)2=E(B).=0, E(B)y= (e/@)(1—8°)-+; 
H(B).=H(B),=0, 
H(B).= —(Be/d")(1—6*)-*, 


in which B=v/c. Thus the components of the 
forces on A and B must be 


F(A)2= F(A).=0, F(A) =eE(A)y—e8H(A)« 
=-(¢/@)(1-B), (5 

F(B)s= F(B).=0, F(B)y=eE(B),—eBH(B). 
= (e8/@)(1— 8"). 


The components of the accelerations of these 
charges are 


a(A)z=a(A),=0, a(A)y= F(A)y(1/m)(1—6*)! 

== (¢/md\(1- 8), 
a(B).=a(B),=0, a(B),= F(B),(1/m)(1 — p)} 

= — (e?/md*)(1—6*), 


where m is the rest mass of each charge, or the 
mass with respect to an observer moving with 
the charges, and m(1—8*)~} is the relativistic 
mass. 

From the point of view of an observer moving 
with the charges, however, the only forces acting 
on the charges are electrostatic. Their compo- 
nents are 


F'(A).=0, F'(A)y= —€/d, F'(A).=0; 
F'(B).=0, F'(B),=e/d?, F’(B).=0. 


Therefore, the components of the accelerations 
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are 


a’(A)z=0, a’(A),= —e?/md?, a’(A),=0, (4°) 
a’(B).=0, a’(B),=e?/md?, a’(B),=0. 


According to the special theory of relativity, 
if the acceleration of a particle with respect to 
one observer is a, and the acceleration of the 
same particle with respect to a second observer, 
in uniform motion with respect to the first, is 
a’, and further, if the particle happens to be 
instantaneously at rest relative to the second 
observer, than @ and a’ are connected by the 
transformation 


a.=az'(1—w/2), ay=a'y(1—w*/e), 
a,=a,/(1—w*/c?). (5) 


Here w is the relative velocity of the two 
observers. 

A comparison of Eqs. (4) and (4’) shows that 
in them the accelerations a and a’ are related in 
exactly the manner stated in Eq. (5). Thus the 
derivation of Eqs. (4) and (4’), based on rela- 
tivistic electrodynamics, leads to a result con- 
sistent with relativistic mechanics and no para- 
dox appears. The point to be emphasized is that 
it is only because relativistic electrodynamics is 
based on action propagated with finite velocity 
c that this agreement is possible. 

The author acknowledges with thanks the 
considerable help he has received from Dr. 
William T. Payne in achieving a straightforward 
and clear presentation of this topic. 


Photosensitive Glass 


KaRLEM RIEss, WALTER C. BoscH, AND T. Topp REBOUL* 
Tulane University, New Orleans, Louisiana 


HOTOSENSITIVE glass is a new type of 
glass in which a permanent photographic 
image of exceptional fidelity is reproduced after 
exposure to ultraviolet light and development 
by heat. The undeveloped glass is clear, of high 
quality, and not distinguishable by its appear- 
ance from ordinary glass. There is no evidence 
of particles in the glass when examined under a 
microscope. The glass may be made in any de- 
sired shape consistent with ordinary commercial 
methods, for example, by drawing, rolling, blow- 
ing, or pressing. The chemical composition of the 
glass is covered by patents. 

Two steps are required to impart a photo- 
graphic image to the glass: first, exposure to 
ultraviolet light through an ordinary film or 
glass negative; second, development at high 
temperature. The light source must produce 
radiation in the 310-340 millimicron region of 
the ultraviolet. After exposure, the glass is still 
clear and transparent. The latent image is de- 
veloped by simply heating the glass uniformly 
to approximately 565°C. Development may be 
completed at any time after exposure provided 


* Present address: Department of Physics, University 
of Pennsylvania, Philadelphia, Pennsylvania. 


the glass is protected from ultraviolet. Times of 
exposure and development can be varied to 
control the depth of penetration and resultant 
color.! 

The method of development of the image, the 
transmission properties of the glass, density and 
exposure time relations were studied. These will 
be considered in the following sections. 


Development 


A circular sample of the photosensitive glass 
was placed on a porcelain crucible, in a hori- 
zontal position, and developed in an electric 
furnace for 30 minutes at 565°C. Temperature 
measurements were made with a _ chromel- 
alumel thermocouple, with the hot junction 5 
mm above the middle of the specimen. As a 
result of this treatment the sample was badly 
distorted. To overcome the difficulty the glass 
was placed on a carbon slab, raised about 3 cm 
above the floor of the furnace, and then de- 
veloped. 

The recommended development temperature 
of 630°C was found to be too high. At this 


1 Anon., ‘‘Photosensitive Glass Developed by Corning,”’ 
Chem. and Eng. News 25, 1822 (1947). 
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Fic. 1. Curves of wavelength versus percentage of 
transmission for unexposed and exposed glass. 


temperature the glass became too soft to hold 
its shape. 


Transmission Properties 


The transmission properties of photosensitive 
glass were studied using tungsten lamp, mer- 
cury arc, carbon arc and hydrogen arc sources. 
Data for the transmission curves shown in Fig. 1 
in the visible region from 400 to 700 millimicrons 
were obtained using a Cenco-Sheard spectro- 
photelometer. The associated galvanometer was 
mounted on a Julius suspension. The source was 
a simple projector with a 6-volt, 18-ampere 
tungsten lamp, power being received from a 
voltage-regulated power supply so that the il- 
lumination could be held constant during a 
series of readings. A converging lens, 7.6 cm in 
diameter, focused the light on the slit. The out- 
put of the photoelectric cell in the spectro- 
photelometer was connected through a shunt 
box to the galvanometer. Because of the pre- 
dominance of yellow and red in the spectrum of 
the lamp, and the high sensitivity of the photo- 
electric cell in that region, it was necessary to 
reduce the voltage applied to the light source 
when using long wavelengths. Transmission 
characteristics of the glass in the ultraviolet 
region were determined by a Beckman spectro- 
photometer, with hydrogen arc source. 

The transmission curves reproduced in Fig. 1 
indicate that after development the glass ex- 
posed the longer time to ultraviolet light is not 
only more dense than the less exposed glass, but 
has an increasing transparency at the longer 
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wavelengths as well as a gentle maximum in the 
400-millimicron region. Two noticeable changes 
were found after exposing and developing. The 
glass which was exposed for a shorter time took 
on a bluish tinge. The glass exposed for a longer 
period became reddish. On closer examination 
by viewing edgewise it was found that the well 
exposed glass did show a blue layer, but this 
layer had penetrated farther into the glass. The 
blue layer was at the edge of the red coloration. 
As the exposure time was increased the thick- 
ness of the red layer in the glass increased. The 
thin blue layer was the same thickness in each 
case, but moved into the glass. 

The thickness of the blue layer may be in- 
creased by exposing the glass to radiation of 
longer wavelength. For example, if an intense 
carbon arc is used and if the very short wave- 
lengths (below 320 millimicrons) are filtered out, 
the red and blue colorations will still be present, 
but the thickness of the blue layer will be 
greater.” 

A sample of the glass was exposed to a mer- 
cury arc-quartz spectrograph combination for 
8 hr and then developed at 565°C for 1000 
sec. Only one broad spectrum line was vis- 
ible, and this was very faint. A longer ex- 
posure time of 23 hr was used to try to 
strengthen this line, and to bring out other lines. 
After development the same line was present, 
with greater intensity this time. No other lines 
were visible, as Fig. 2 shows. The line affecting 
the glass was in the 312-313 millimicron region 





Fic. 2. Sample of photosensitive glass showing the 312-313 
millimicron line of the mercury arc spectrum. 


2R. H. Mariner, private correspondence from Corning 
Glass Works. 
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of the arc spectrum. The 312 and 313 milli- 
micron lines could not be resolved with the 
instrument used. 

The sensitivity curve of the glass has a peak 
at about 320 millimicrons.? On either side the 
302- and the 334-millimicron lines are already 
too far away from the sensitivity peak to affect 
the glass. Since this spacing of strong lines is 
characteristic of the mercury arc, there can be 
no control of the blue layer thickness with this 
source. It is, however, evident that the region of 
very low transmission and the region of sensi- 
tivity of the glass coincide. 

When a small carbon arc and quartz spectro- 
graph were used, and the glass was exposed for 
22 hr, no results were observed when the glass 
was developed. 

The curves in Fig. 1 show that the percentage 
transmission for the three glass samples falls 
off rapidly at short wavelengths, beginning close 
to the 360-millimicron region. At 310 milli- 
microns transmission is barely perceptible even 
with the two least exposed glasses. The well 


exposed glass absorbs all wavelengths below 340 
millimicrons. 


An experiment was done with a hydrogen arc 
whose spectrum extended to shorter wave- 
lengths than these limits. No light was trans- 
mitted through the glass in the range from 310 
millimicrons to 200 millimicrons. 


Influence of Development Time 
and Temperature 


The photosensitive glass was exposed for 40 
min to the mercury arc at a distance of 45 cm. 
This time was sufficient for satisfactory contrast, 
but not excessive, and was chosen to avoid any 
variations which might arise as a result of a 
reversal effect. 

White light from the projector was sent 
through a No. 244 Corning filter and then 
through the photosensitive glass while it was 
developing in the furnace. The emergent light 
was received on a photovoltaic cell, permitting 
an observation of the change in transmission 
characteristics during developing. Figure 3 shows 
the results in graphical form. 

There was no change in the transmission char- 
acteristics of the glass for the first 3500 sec, 
while the furnace was warming up. The glass 
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Fic. 3. Curves showing the relationship between trans- 
mission and development time, and between furnace tem- 
perature and development time. 


showed no signs of developing until the tempera- 
ture had reached 540°C. Since the prescribed 
development temperature was far higher than 
this, an arbitrary chosen constant furnace tem- 
perature of 565°C was maintained. The tempera- 
ture was held constant until there was no further 
change in the transmission characteristic of the 
glass. The total time of development after the 
initial warming up period was 2750 sec. After 
development the temperature was lowered, but 
the transmission characteristic remained con- 
stant. ; 


Variation of Density with Time of Exposure and 
Time of Development 


Tests were made to see how different de- 
velopment times would change the density of 
variously exposed pieces of glass. The time of 
development was taken to be that time during 
which the temperature in the furnace was be- 
tween 550 and 565°C. The glass samples were 
exposed at a distance of 45 cm from the mercury 
arc lamp for 5, 10, 20, 40, 80 and 160 min. 
After exposure the glass samples were placed in 
a special holder of brass tubing, separated by 
metal rings, and held secure by a large over- 
sized ring. This made possible the simultaneous 
development of several samples under nearly 
identical conditions. The specimens were heated 
for a prearranged time, removed from the fur- 
nace and the density of each glass determined. 
The ends of the furnace were sealed to eliminate 
circulation of air during the development. The 
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Fic. 4. Density as a function of exposure for 
several development times. 








samples were then replaced in the tube and de- 
veloped again. This was repeated five times, so 
that there were six development times for a 
total of 2600 sec. 


The relation used to determine density was 
D=logiolo/T, 


where J is the intensity of light transmitted 
through an unexposed piece of glass, and J the 
intensity of light transmitted through the ex- 
posed and developed glass. The intensity read- 
ings represented in Fig. 4 were made using a 
photovoltaic cell with a broad beam of emergent 
light covering the entire face of the cell. 















































































































Fic. 5. Sample of photosensitive glass containing an 
image from a lantern slide negative. 
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One of the important relations characteristic 
of photosensitive materials is that between the 
time of exposure and the density. The normal 
curve for a photographic film may be found in 
standard text books on photography as a plot 
of density against logiot, where ¢ is the exposure 
time.* The curve obtained with photosensitive 
glass may be compared with that of a normal 
photosensitive emulsion. Figure 4 gives curves 
characteristic of various stages of development. 
Curves A, B and C do not follow the normal 
curve, for the density continues to increase in 
the region of abnormal over-exposure instead of 
decreasing. Curve F represents the totally de- 
veloped glass. 

For a photographic emulsion a low value of 
density occurs upon development when the ma- 
terial has not been exposed. This is called fog. 
A sample of the photosensitive glass which had 
not been exposed was heated for 3000 sec at a 
temperature of 565-600°C. After cooling there 
were no signs of fog, since in this case the 
density remains zero after development. The 
underexposed region for the glass extends from 
about ¢=5 to t=10 min. This portion of the 
curve is similar to a normal curve in the under 
exposed region. The region of correct exposure is 
that portion of the curve between ¢=10 and 
t=15 min. On the logarithmic scale of Fig. 4 it 
appears rather narrow. 

From ¢=15 to t=80 min is the region of over- 
exposure. If compared with the curve of a normal 
photosensitive emulsion, it is seen that this region 
covers a much greater range on the exposure 
time coordinate than is the case with the normal 
photographic emulsion. This difference may be 
due to the fact that the exposed and developed 
glass has two layers of color, red and blue, as 
mentioned above. Initially, with small exposures, 
the curve is normal since only the blue color is 
present. However, with longer exposures the red 
color begins to come into the glass and the com- 
bination of these two colors gives rise to the 
broad region. 

The region beyond ¢=80 min is the abnor- 
mally overexposed region, where partially de- 

3C. E. Kenneth Mees, Theory of the photographic process 
(Macmillan, 1944); Keith Henney and Beverly Dudley, 
Handbook of photography (Whittlesey House, 1939); J. E. 


Mack and M. J. Martin, The photographic process 
(McGraw-Hill, 1939). 
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veloped samples show an increase in density, 
and a totally developed sample, as shown by 
curve F, Fig. 4, shows a decrease in density as a 
result of very long exposures. This reversal 
effect is the same as that given by a normal 
photographic emulsion. The time of develop- 
ment was chosen so that there would be a marked 
change between successive curves in Fig. 4. 
Beyond 2000 sec the effects of further develop- 
ment were hardly noticeable. 

The gradient y, or slope of the curve in the 
correct exposure region, is an important defining 
characteristic of a photosensitive material. 


y=dD/d (logit), 


where dD is the incremental change of density 
with the corresponding incremental change of 


logiot. Values of y for the glass samples are listed 
in Table I. 


Photographic Prints 


The mechanism of image production is not 
clearly understood. Apparently, submicroscopic 
metallic particles are present in solution in the 
glass mix. The incident ultraviolet light pre- 
cipitates these particles, and they become 
colored. A three-dimensional effect is observed 
because the shadowed areas penetrate farther 
into the glass than the highlighted areas. Re- 
production of extremely fine detail may be due 
to the lack of grain in the glass, and to the small- 
ness of the particles making up the image. The 


TABLE I. Values of y for photosensitive glass. 


Time of develop- 
ment (sec) 


200 400 600 1000 
Gamma 


1600 2600 
0.44 0.80 0.85 0.98 


1.41 2.40 


metallic particles are much smaller than the 
silver grains in photographic emulsions.! 

An ordinary Ansco photographic film negative 
was placed over the glass and the combination 
exposed for 10 hr to a mercury arc at a distance 
of 45 cm. When the glass was developed it was 
found that it had not been affected by this ex- 
posure. The film had absorbed all of the effective 
radiation. Other photographic films do not have 
this absorption property to nearly so great an 
extent. 

It was well known that lantern slide glass 
plates transmit the ultraviolet light which affects 
the photosensitive glass. A lantern slide negative 
was, therefore, bound to the photosensitive 
glass and exposed for 15 hr. The glass was de- 
veloped into a clear, distinct picture with a 
three dimensional effect, as shown in Fig. 5. The 
contrast in this reproduction is due to the red 
color of the exposed glass and to the use of 
orthochromatic film in copying it for publication. 
No filters were used. 

The samples of photosensitive glass used in 
this investigation were obtained through the 
courtesy of the Corning Glass Works. Investiga- 
tions of other optical properties of the glass are 
in progress. 


One thing that never seems to occur to some of those who worry so loudly about the discretion of 


scientists is that the information in question, in many instances, was not given by the government * 
to the scientist in the first place. Rather it was given by the scientist to his government. And as 
many of us know from personal experience it was sometimes most difficult to get the government to 
listen. 

The history of the atomic bomb project is an interesting case in point. News of the discovery of 
uranium fission reached this country in January, 1939 from Germany. Within a short time quite 
a few American physicists recognized the possibility of useful release of atomic energy and of 
making an atomic bomb. Then started a process of trying to interest the government with no apparent 
action resulting. We physicists in the meantime voluntarily adopted secrecy policies which kept 
this information from the public and from other countries. 

After some months of frustration a direct appeal to the President was made and he saw to it that a 
program of work was started under the general supervision of my predecessor, Dr. L. J. Briggs. 
The secret was so well kept that most of the staff of the National Bureau of Standards were unaware 


of the existence of an atomic bomb project before the official announcement was made,—E. U, 
Conpbon (1948), 





Ancient Science in the Modern Curriculum 


J. J. G. McCue 
Smith College, Northampton, Massachusetts 


HE present quickened interest in physics as 

a cultural subject gives some teachers of 
physics not only an opportunity, but an obliga- 
tion, to consider with their classes the nature and 
history of scientific activity. Much light is 
thrown on these matters by a consideration of 
Greek and Roman science. A very little reading 
in the subject is enough to convince one that 
Greek science was far from being the contemptible 
guesswork that it is often supposed to have been. 
By inquiring into the science of the classical 
period we can see things about science which 
might not be obvious from a study of the modern 
world alone, but which are of the utmost im- 
portance in the present period of uncertainty 
concerning the relation of science to society. The 
purpose of this article is to outline a few of the 
accomplishments of the Greeks, showing the 
extent to which physical science developed in 
their hands, and to point out some of the 
inferences that can be drawn from the rise and 
fall of science in the ancient world. 

We may begin by observing that the Greeks 
invented science. At least until the advent of 
quantum mechanics, science has rested on the 
hypothesis of ‘‘a constant and universal sequence 
of cause and effect.’’ The Greeks were the first to 
believe in this hypothesis. To be sure, not all 
Greeks believed in it, but not all Americans be- 
lieve in it today, many of our religious practices 
being at variance with it. The important thing is 
that an appreciable number of Greeks believed in 
the hypothesis, while nobody in earlier cultures 
seems to have done so. 

As far as method is concerned, the fundamental 
characteristics of physical science are (i) obser- 
vation or experimentation and (ii) expression of 
findings in mathematical form. We take these 
tremendous intellectual inventions for granted 
now, but if the Greeks had bequeathed nothing 
else to our civilization, we should still owe them 
eternal respect for this legacy. One may ask, 
however, whether they did anything with science 
after they had invented the principles on which 
it rests, and whether we are not justified in re- 


garding them as mere speculators rather than as 
scientists. Before summarizing a few of the actual 
accomplishments of the Greek scientists, it will 
be well to mention some of the reasons why Greek 
science seems meager to us. 

The first and most important reason is that, as 
Sarton! puts it, science and technology are almost 
the only human activities which are cumulative. 
Failure to recognize this difference between 
science, in which the work of many generations is 
at our disposal, and the arts, in which each gener- 
ation fends for itself, leads us to underestimate 
Greek accomplishment in science. Second, our 
impression of Greek science is too often based on 
the works of Plato and Aristotle, whose chief 
interests and skills lay in other fields of learning 
than physical science. Third, the literary tradi- 
tion has had such a hold on classical scholars that 
they have almost completely ignored the Greek 
scientists. Finally, the works of the best Greek 
scientists were at variance with the religious 
views of the Christians, Jews and Moslems on 
whom we are dependent for the transmission of 
Greek learning. Hence our knowledge of these 
works is fragmentary. 

Let us now consider some of what the ancient 
scientists actually accomplished, casting an oc- 
casional glance at the relationship between 
science and society during this period.? To review 
all of the scientific work of the Greeks is of course 
impossible and unnecessary for the present pur- 
pose. By considering only the progress and decline 
of the ancient world’s knowledge of the solar 
system, including the earth, we shall see that the 
Greeks, using observation, experiment, and math- 
ematics, gained a very considerable store of posi- 
tive knowledge, that the Greek scientific enter- 
prise succeeded to a degree that commands 
respect even in our time, and that it then 
withered away. By examining the causes of the 


1G. Sarton, Introduction to the history of science (Carnegie 
Institution of Washington, 1927), vol. I, p. 4. 

2 The principal references used in preparing this summary 
have been Sarton (reference 1) and G. Singer, A short 
history of science (Oxford, 1941). Singer’s excellent book can 


be recommended to anyone seeking an orientation in the 
subject. 
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withering, we may hope to throw light on prob- 
lems of today. 

Astronomical records, used for predicting 
the proper time to plant crops, existed in 
Mesopotamia at a very early date. The annual 
flooding of the Nile valley is said to have given 
rise to the science, or perhaps more properly the 
technology, of surveying; at any rate, “‘geome- 
try’ is Greek for ‘‘earth-measurement.” These 
pre-Greek efforts yielded a considerable accumu- 
lation of recorded data and rules of thumb, but 
science as an effort to correlate phenomena by 
means of universally valid principles seems cer- 
tainly to have been invented by the Greeks. In 
fact, its invention is usually attributed to a par- 
ticular Greek, Thales of Miletus,* who flourished 
about 580 B.C. Thales discovered some geo- 
metrical generalizations, such as some of the 
properties of similar triangles. These properties 
he used to solve practical problems, for example, 
finding the distance from shore of a ship at sea. 

Thales also thought about the nature of things, 
reaching the conclusion that everything is made 
of water. This is naive, but not silly; it was not 
long ago that we thought everything was made of 
hydrogen. The important point is that Thales 
tried to explain the universe in rational terms. 

Pythagoras (530 B.C.)* usually considered the 
inventor of mathematical proofs, also founded 
mathematical physics with his discovery of the 
simple ratios of the lengths of strings sounding 
the musical intervals (when the strings all have 
the same tension and linear density). This dis- 
covery must have been based on experiment. With 
mathematical proofs, the idea of experimenting, 
and the idea of applying mathematics to the 
results of experiments, science in the modern 
sense was fairly started. 


While the Pythagorean school was flowering in’ 


southern Italy, the Miletan school’ founded by 


3 An ancient city on the west coast of what is now 
Turkey. 

4 The dates of important men are mentioned in order to 
establish a rough scale of time. They indicate the year when 
the man was about forty. 

5 The reader should perhaps be warned at the outset 
against one of the primary misconceptions about Greek 
science: that it was a sporadic activity, engaged in from 
time to time by men working in isolation. It will become 
clear that Greek science was actually a continuous activity, 
carried on in schools that endured for generations by men 
who were cognizant of the activities of other workers. The 
misconception doubtless arises from the random and 
fragmentary nature of the documents which happen to 


Thales was still in action. One of Thales’ dis- 
ciples, Anaximander (570 B.C.), used a sun dial 
to study the motions of the sun, and to determine 
the solstices and equinoxes. Another, Anaximenes 
(530 B.C.), reached the conclusion that the moon 
shines by reflecting light that comes from the sun. 

About 450 B.C. wealthy Miletus, the greatest 
Greek city of its time, was conquered by the 
Persians, and the intellectual center of Greece 
shifted from Asia Minor to Athens. Anaxagoras, 
a contemporary and friend of Pericles and 
Euripedes in Athens, thought the earth to be flat, 
but correctly explained eclipses of the sun and the 
moon. For teaching that the sun is a huge ball of 
white-hot stone® he was charged with impiety and 
is said to have found it wise to return to Asia 
Minor. 

Anaxagoras’ contemporary, Socrates (430 
B.C.), must be mentioned because of his two- 
fold influence on the progress of science, which he 
aided by combating loose reasoning, but impeded 
by turning men’s attention to ethical problems, 
with which he himself was preoccupied. 

Socrates’ most distinguished pupil, Plato (390 
B.C.), was also concerned largely with problems 
of conduct, but contributed much to the logical 
structure of mathematics and fostered the study 
of astronomy. His disciple, Eudoxus (370 B.C.), 
put astronomy for the first time on a mathe- 
matical basis, developing a very ingenious system 
of concentric spheres to account for the retro, 
gressions in the motions of the planets. It is im- 
portant to mention that he did not regard these 
spheres as material. He also measured the length 
of the solar year, making an error of only about 
one minute. 

The idea that the earth is a sphere dates back 
to Pythagoras; it seems to have won general ac- 
ceptance at least as early as the generation after 
Plato. Heraclides, another of Plato’s pupils, 
postulated the rotation of the earth to account 
for the apparent diurnal rotation of the heavens, 
and discovered that Venus and Mercury—the 
two planets whose orbits are smaller than that of 
the earth—revolve about the sun. 

Aristotle (340 B.C.), also a pupil of Plato, 


have survived the long period when they were regarded 
without interest or with disfavor. 

6 He seems to have based this judgment on a huge 
meteorite which fell near the Hellespont. See Sarton, refer- 
ence 1, p. 86. 
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understood the function of research and con- 
tributed enormously to biology by careful in- 
vestigation, but his work in physics was less 
fortunate. From our point of view, he made two 
major mistakes: he attempted a synthesis ex- 
plaining all dynamical phenomena in terms of 
first causes, without having a sufficiently detailed 
knowledge of the phenomena, and he underesti- 
mated the importance of quantitative observa- 
tion. His extensive writings on physics were very 
influential during the Revival of Learning; they 
contribute heavily to the Greeks’ reputation as 
poor scientists. 

After the almost coincident deaths of Aristotle 
and Alexander the Great, Alexandria took the 
place of Athens as the scientific center of the 
western world. More significant than the change 
in site was the accompanying change in method. 
Scientists became a specialized group, no longer 
seeking to embrace all knowledge and study the 
universe as a whole. 

Perhaps the greatest scientist at Alexandria 
was Aristarchus of Samos (270 B.C.), who ranged 
through all of the mathematical sciences of his 
time, but is chiefly noted for his works on as- 
tronomy, of which the only survivor is On the 
sizes and distances of the sun and the moon.’ The 
relative distances of the sun and the moon from 
the earth were determined by estimating the 
angle A between lines of sight to the sun and the 
moon when the moon was exactly at the half 
(Fig. 1). Aristarchus took A to be 87°, and com- 
puted the ratio of the desired distances. Since A 
is actually about 89° 50’, Aristarchus’ result was 


MOON 
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Fic. 1. Illustrating Aristarchus’ method of finding the 
relative distances of the sun and the moon. 


7 For an English translation of this work, see Sir Thomas 
Heath, Aristarchus of Samos (Oxford, 1913), which also 
contains a summary of Greek astronomical progress before 
Aristarchus. 
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much too small—about 19 instead of 346. He 









the 
knew from watching an eclipse that the sun and a 
the moon subtend very nearly the same angle at § bette 
the earth, and therefore was able to reason that § of th 
the diameter of the sun must be about 19 times J phen 
that of the moon, which made the volume of the § elabo 
sun about 7000 times that of the moon. More- § centr 
over, estimating the angle subtended by the § tical 
moon at the earth to be® 2° and the diameter of § impr 
the earth’s shadow on the moon during a lunar J mind 
eclipse to be twice the diameter of the moon, he § caus¢ 
deduced theorems to show the volumes of the J weig’ 





moon, earth and sun to be as 0.05, 1 and 300. 
These figures are simplifications of those of 
Aristarchus, who, lacking trigonometric tables, 
calculated lower and upper bounds of the ratios. 

Because of low accuracy in the measurements, 
these results are far from correct, but they are of 
great intellectual significance because they 
showed that while the moon is smaller than the 
earth, the sun is very much larger. It was perhaps 
this investigation that led Aristarchus to his 
crowning achievement, which was nothing less 
than the detailed anticipation of the Copernican 
hypothesis, in a book which is now lost. Several 
predecessors of Aristarchus had taught that the 
earth rotates on its axis, and Heraclides, seventy 
years or so earlier, had found that Mercury and 
Venus revolve about the sun. We are told by 
Archimedes, a younger contemporary, that 
Aristarchus took the sun to be one of the fixed 
stars, about which the earth and the other planets 
revolve in circular orbits; the moon revolves 
about the earth. He stated that the diameter of 
the earth’s orbit is negligible in comparison with 
the distances to the stars. Therefore no parallax 
was to be expected, and the major objection to 
the motion of the earth (raised by Aristotle*) was 
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Aristarchus’ system and mentioned it in a pas- 
sage which he subsequently deleted.!° 

The successors of Aristarchus improved on his 
measurements of the sizes and distances of the 
moon and the sun, but his heliocentric theory was 
soon rejected, probably because it did not ac- 
count accurately for the motions of the moon and 


kn 
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8In a presumably later work, Aristarchus announced 
that this angle is 30’, a much more accurate value. 

9 De coels ii, 14 

10 Heath, reference 7, p. 301. 
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the planets when improved observations, made 
with improved instruments, became available. As 
better observations brought to light more details 
of these motions, the astronomers “saved the 
phenomena,” or made the theory fit the facts, by 
elaborating the number of epicycles in the geo- 
centric system. Not until Kepler suggested ellip- 
tical orbits was the heliocentric theory a great 
improvement on the geocentric one. One is re- 
minded of the rejection of Prout’s hypothesis be- 
cause of improved measurements of atomic 
weights, and its resurgence after the discovery of 
isotopes. 

To go from Thales’ idea of the earth floating on 
the water like a cork to Aristarchus’ heliocentric 
theory took the Greeks just about three centuries. 
Among the astronomers who followed Aristarchus 
may be mentioned Hipparchus (150 B.C.), who is 
notable for his invention and improvement of 
astronomical instruments and his discoveries of a 
nova and the precession of the equinoxes, and 
Ptolemy (150 A.D.), whose chief accomplishment 
was the collection and exposition of Hellenistic 
knowledge of the heavens and the earth in 
his Almagest and Geography. The quality of 
Hipparchus’ work may be judged from his value" 
of the mean lunar month, which differs from the 
present value by less than 1 sec. Ptolemy’s works 
have been passed on to us by the Arabs, but were 
not translated into Latin until the close of the 
Middle Ages. Hipparchus worked partly, and 
Ptolemy wholly, at Alexandria ; some of Ptolemy’s 
work was original, but he depended very largely 
on his predecessors, particularly Hipparchus. In 
spite of living three centuries apart, these two 
men were almost collaborators. We therefore 
notice a striking decrease in the rate of progress, 
in comparison with the earlier period. 

It is instructive to consider also the progress of 
knowledge concerning the earth itself. The work 
of Eratosthenes (235 B.C., a generation after 
Aristarchus) lay primarily in this field. His meas- 
urement of the diameter of the earth is an ex- 
ample of Greek science at its best. At Syene 
(Assuan) in Egypt, at noon on the summer 
solstice, a vertical post cast no shadow, indicating 
that the sun was directly overhead. Eratosthenes 
measured the length of the shadow cast at the 


| Cambridge ancient history (Cambridge, 1928), vol. 7, p. 
311, - 
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Fic. 2. Illustrating Eratosthenes’ method of finding the 
diameter of the earth. 


same time by a vertical rod at Alexandria, which 
he took to be 5000 stadia north of Syene. The 
lengths of the rod and its shadow gave the angle 
6 (Fig. 2) between the earth’s radii to the two 
points, thereby enabling him to compute the 
diameter of the earth.!2 His measurement of 6 
was very good, but his value for the base line was 
1 percent too small, and his value of 7820 statute 
miles! for the diameter of the earth is therefore 
1 percent smaller than the true distance from 
pole to pole. 

Eratosthenes also made a map of the world as 
he knew it, and suggested that only the width of 
the ocean prevented the circumnavigation of the 
globe. He estimated the breadth of Europe and 
Asia to be 7600 miles; actually the greatest dis- 
tance from the European coast to the Asian 
coast, along a parallel of latitude, is about 7640 
miles (Lisbon-Fusan).'* He took Ceylon and the 
Shetland Islands to be the southern and northern 
limits of the inhabited world, and placed the 
Shetlands 3700 miles farther north than Geylon. 
The error in this distance cannot be computed 
without knowing which part of the Shetlands he 


12 For the practical details of this experiment, see Pauly- 
Wissowa, Real-Encyclopadie der Classischen Altertums- 
wissenschaft (Stuttgart, 1909), vol. 6, pp. 363-366. 

18 The length of his stadium is not known with certainty, 
but is thought to be 0.0977 statute miles (see Heath, refer- 
ence 7, p. 339). The startling accuracy of all three of the 
measurements quoted here gives some confidence in the 
correctness of this conversion factor. 

14 [t is possible, but unlikely, that Eratosthenes had in 
mind the longest great-circle arc between the coasts. The 
true value for this distance is 7800 miles (Lisbon-Singapore). 
In any case, if Columbus had known this distance and the 
diameter of the globe as well as Eratosthenes did, he would 
scarcely have set out to cross the ocean. 
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referred to, but the southern limit of Ceylon is 
3750 miles south of the center of Mainland, the 
chief Shetland isle. 

Four centuries later, Ptolemy’s Geography 
showed some improvement in geographical know]l- 
edge, but here too the rate of progress had obvi- 
ously fallen off. We may say that Ptolemy’s map 
of the world, reconstructed from his latitudes and 
longitudes, represents correctly the broad fea- 
tures of northern Africa, southern Asia as far as 
Indo-China, and Europe exclusive of Scandi- 
navia. In the roughly contemporary Roman 
world, however, Tacitus, writing on geographical 
matters, described the earth as flat, and by 600 
A.D., when Isidore of Seville compiled his 
encyclopedia, the map of the world as far as 
Europeans were concerned had degenerated so 
much as to be scarcely recognizable. 

We see, therefore, that from 550 to about 150 
B.C., the Greek world developed to a high degree 
the technic of scientific investigation of certain 
types of problems. We have mentioned cosmog- 
raphy and geography ; very considerable progress 
was also made in some branches of biological 
science, where dissection of animals and men led 
to anatomical discoveries which were in some 
cases not repeated until the nineteenth century. 
We naturally inquire why this scientific activity 
died out. The complete answer is of course very 
complex, but a few of the reasons may be 
mentioned. 

When the empire of Alexander fell apart after 
his death, the unsettled political conditions were 
unfavorable to the development of science. The 
lot of most people was so miserable that philoso- 
‘phy concerned itself primarily with questions of 
ethics and religion; science seemed unimportant. 
In Egypt, however, stability and wealth per- 
mitted the ruling Greeks to foster the develop- 
ment of science by the foundation and support of 
an institute like a university, the Museum at 
Alexandria. Scholars and scientists came here 
from all over the Greek world, and ancient 
science entered on its golden age. These scientists 
were specialists, not so specialized as the modern 
scientist, but much more so than Plato and 
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Aristotle. In the second century B.C. Rome 
established control over the Eastern Mediter- 
ranean, although Egypt retained nominal inde- 
pendence until about 50 B.C. The educated 
Romans tried to assimilate Greek culture, and 
succeeded to some extent in philosophy, literature 
and the fine arts, but they failed completely in 
science. Among the reasons for this were un- 
doubtedly the following. 

(i) Science failed to appeal strongly to them 
because of their religious convictions, which em- 
phasized the importance of conduct and mini- 
mized the importance of understanding the ma- 
terial world. Partly, perhaps largely, for military 
reasons, they did have a strong interest in 
technology, but this of course is not science. The 
validity of the distinction is confirmed by the 
difference in results. 

(ii) They never mastered the higher branches 
of Greek mathematics, and therefore could not 
produce original work in the exact sciences. This 
is probably a consequence of (i). 

(iii) Above all, they failed to recognize that 
science is an operation, rather than a collection of 
facts. They tried to take over and preserve the 
results of Greek science. This attempt failed be- 
cause the body of real knowledge inevitably be- 
came more garbled and diluted with each genera- 
tion, so that it was almost entirely lost long before 
the fall of Rome. 

The fact that science is an operation, rather 
than a collection of facts, is probably the most 
important thing the teacher of physics as a 
cultural subject has to impress on his present 
classes, if their preconceived notions of science 
are to be replaced by real understanding. Starting 
the elementary physics course with a lecture or 
two on ancient science affords a convincing 
demonstration, based on simple subject matter, 
that science is not a thing that inevitably pro- 
gresses. It does not grow like a weed. On the 
contrary, it is a delicate organism depending in a 
very sensitive way on its environment. !t has 
died before, and without the proper conditions of 
life it will die again. 
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Principal Series of Sodium in Absorption 


J. RAND MCNALLY, Jr.* 
Massachusetts Institute of Technology, Cambridge 39, Massachusetts 


HIS note describes the results of a qualitative in- 
vestigation into the principal series absorption’ in 
sodium vapor, using an absorption tube designed pri- 
marily for use by students in the third-year course in 
physics at the Massachusetts Institute of Technology. 
The apparatus is described elsewhere,! but essentially 
consists of a ten-inch seamless stainless steel tube, sur- 
rounded by a heater coil and insulating material, and hav- 
ing water-cooled ends to which quartz windows are 


attached. Stockbarger’s article describes the tube’s possi- 
bilities as follows. 


Apparently the number of principal series members 
of sodium which can be observed with the apparatus is 
determined largely by the resolution of the optical sys- 
tem. With a mercury arc source, students obtain about 
15 members with a small quartz spectrograph and 30 
or more with a medium size instrument, - - - 


In the present work clean sodium metal was placed in 
the apparatus which was then outgassed at a high tempera- 
ture. A 35-foot grating spectrograph, having a dispersion 
of 0.76A/mm in the first order and a resolving power ex- 
ceeding 100,000, was used to obtain absorption by the 
sodium at temperatures above 600°C out of the high in- 
tensity continuum in the spec- 
trum of a quartz pool-type 
mercury arc. The series from 
3S—15P(2444.195A) to the 
limit at about 2410A is repro- 
duced in Fig. 1. The last series 
member detected on the origi- 
nal plates was the 3S—52P 
transition, which indicates a 
resolving limit of approxi- 
mately 0.095A. Resonance 
broadening and broadening due 
to residual gases are probably 
the factors that limit the num- 
ber of series members attainable. 

The Stockbarger absorption 
tube is thus capable of resolving 
the principal series in sodium 
to a degree which compares 
very favorably with the early 
work of Wood.? Until very re- 
cently,*# Wood's observation of 
the 3S—59P electron transition 
involved the highest quantum 
state ever observed. It is of in- 
terest to point out that this was 


—2444.195A 


Fic, 1. High members of the principal 
series of sodium in absorption. 


observed with a 2.8-meter absorption tube, and that the 
limit of observation was probably set by the resolving 
power of his spectrograph since the last twelve lines were 
observed in an interval of less than one angstrom. 


* Now at Carbide and Carbon Chemical Corporation, Isotopes 
Branch, Oak Ridge, Tennessee. 

1D. C. Stockbarger, J. Opt. Soc. Am. 30, 362 (1940). 

2R. W. Wood, Astrophys. J. 43, 73 (1916). 

3H. R. Kratz and J. E. Mack, Rev. Mod. Physics 14, 104 (1942). 

4J. R. McNally, Jr., J. P. Molnar, W. J. Hitchcock and N. F. Oliver, 
to be published in J. Opt. Soc. Am. 


On the Area-Moment Propositions of Mechanics 
A. W. SIMON 
University of Tulsa, Tulsa, Oklahoma 


HE proofs of the area-moment propositions in me- 
chanics can be simplified over those usually given 
in the literature, and the propositions themselves can be 
put in a form of more general application, as follows. 
Proof of the first area-moment proposition.—The dif- 
ferential equation of the elastic curve is 


d?y/dx? = M/EI, (1) 


where M is the bending moment, J the moment of inertia 
of the section, and E the modulus of elasticity. Integrating 
between the limits x =a and x=), we have 


bdty b MM 
J. es 


whence, for beams of constant cross section, 


()_.-@)_- 
dx/,.» \dx/,.. EI 


This is equivalent to 
tan@,—tan@,=A/EI. 


Finally, in view of the fact that the angles involved are 
very small, the last equation can be written in the form 


6,—0.=A/EI, (2) 


where, as usual, A is the area between the limits x=a and 
x=b under the corresponding moment curve. 

Proof of the second area-moment proposition.— Multi- 
plying both sides of Eq. (1) by (b—x), and integrating 
between the same limits as above, we have 


f (— ~x)ae= js eel = 


Integrating the left-hand member by parts gives 
Ab*/EI, 


where b* is the distance of the centroid of the area from 
the line x=. Rearranging and taking account of the fact 
that @. is small yields, finally, 


yo=Yat(b—a)0.+Ab*/EI. (3) 


Yo —Ya—(b—a) tande= 
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The last equation reduces to the usual form provided 
¥a=0 and 6,=0. The geometric interpretation of the three 
terms on the right-hand side is obvious. 

It will be noted that the generalized form of the second 
proposition allows the deflection to be obtained in every 
case by purely analytic methods, without resorting to 
additional geometric constructions. In particular, by the 
application of both propositions, a sufficient number of 
algebraic equations can be written (by suitable choice of 
pairs of points a, b) to solve any problem completely. 

Application of the generalized area-moment proposi- 
tions.—Consider the case of any simply supported beam. 
In order to deduce the deflection y at a section located a 
distance x from the left-hand support (x =0), we can apply 
the second proposition twice, once setting a=0, b=x, and 
again setting a=0, b=/; whence we have 


= 0+x00.+A ozdoz*/EI, 
0=0-+16 +A oiboi*/EL. 


Elimination of 69 gives the deflection y directly. 

Alternate form of the second area-moment proposition.— 
If we multiply both sides of Eq. (1) by (x—a) and carry 
out the integration as before, there results 


Yo= Yat (b —a)%+Aa*/EI. 
In view of Eq. (2), this can be written as 
¥o=Yat(b—a)(A/EI+6.)+Aa*/EI, 


where a* now represents the distance of the centroid of 
the area from the left-hand vertical, that is, from the line 
x=a. 


A Dynamic Atom Model 


Joun B. UNDERWOOD 
High School, Grass Valley, California 


HIS Bohr atom model is comparable to a lecture table 
planetarium, in that atomic particles are made to 
describe circular orbits like planets. Obviously no atomic 
model can duplicate proportional dimensions. This is not 
even attempted in an astronomical planetarium. The pur- 
pose of this model is to demonstrate the effect of particle 
motion and to show how it may determine the properties 
of an atom. The nucleus is made up of particles represent- 
ing protons and neutrons, caused to move about in the 
smallest space mechanically possible. When in motion, the 
model gives the impression of a single nuclear mass. About 
the nucleus revolve the planetary electrons in appropriate 
orbits, giving the impression of cloud layers rather than 
particles. Other interesting phenomena such as nuclear 
magnetic effects can be illustrated. But it is our purpose 
here merely to describe the mechanical action of the model. 
The oxygen atom proved ideal for our model. The eight 
protons and the eight neutrons in the nucleus just fill the 
four spokes on each of the four spindles used. The electrons 
are secured at the ends of the wire spokes common to the 
protons, thus illustrating the electrostatic bond of attrac- 
tion between them. The two K orbit electron spokes are 


Fic. 1. The oxygen atom model at rest. 


shorter than those carrying the outer orbit electrons, 
illustrating energy levels. 

The general appearance of the model at rest is shown 
in Fig. 1. A vertical pipe, about four feet long and 
mounted on a tripod base, serves as a supporting stand 
and drive shaft housing. At the top of the pipe is the driv- 
ing head upon which are mounted four spindles rotated by 
miter gears, one of which is driven by a shaft coming up 
through the supporting pipe from the driving gear operated 
by a crank at the base. Attached to each of the four spindles 
are wire spokes that carry the neutrons, protons and elec- 
trons. Hollow rubber balls are used for neutrons and pro- 
tons, and spherical metal buttons for electrons. The white 
balls represent protons and the dark ones neutrons. These 
are free to rotate (see Fig. 2) on the wire spoke axis, thus 
illustrating magnetic effects within the nucleus. The eight 
wire spokes not carrying electrons terminate with the 
neutrons. 

The mechanical details of the head are shown in Fig. 3. 
Standard }-inch pipe, a cross tee, three short nipples, a 
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Fic. 2. The oxygen atom model in motion. 


coupling, and two pipe caps were used for supports as 
shown. It will be seen that the miter gear at the top is 
the driver and that the other gears rotate in unison with 
it. The miter gears were cut from sheet metal and soldered 
on short brass tubing hubs which were made to fit the 
pipe nipples. 

Figure 2 shows the appearance of the atom head when 
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Fic. 3. Details of the driving mechanism of the oxygen atom model. 


the crank at the base is turned two times per second. The 
blur at this speed is mostly photographic, the nuclear 
particles appearing to the eye rather to mill around as we 
would like to show them. At about double this speed the 
blur becomes real to the observer. The electrons appear 
only as circular orbits even at slow speeds. One of the hori- 
zontal orbits in Fig. 2 shows the paths of the two inner 
electrons. The remaining six orbits represent the outer 
electron orbits. In operation the electron orbits stand out 
much more clearly than it is possible to photograph, and 
the apparatus never fails to fascinate spectators. 


RECENT MEETINGS 


Southern California Section 


The Southern California Section of the American Asso- 
ciation of Physics Teachers presented its fourth annual 
competitive physics test for high school students on 
Saturday, May 15, 1948. Examination centers for this 
test were established at the University of California at 
Los Angeles, the University of California at Santa Barbara, 
the University of Redlands, and San Diego State College. 
The examination was two hours in length and included 
definitions of basic quantities and principles, plus problem 
work in the fields of mechanics, properties of matter, heat, 
sound, electricity, and light. The problems were graded on 
the basis of method as well as on accuracy. 

Students participating were nominated for the test by 
their high school physics teachers. There was no limit to 
the number of students any one teacher could nominate. 
This year 193 students representing 47 high schools in the 
Southern California area competed. The fifteen top rank- 
ing students on the test were presented with certificates of 


award from the Section. In addition, duplicate certificates. 


were presented to the schools represented by the winning 
students. These fifteen winners were: 


Jon Mathews, Flintridge Prep. School, Pasadena; Francis Scott, 
John Marshall High, Los Angeles; Bill Colbert, Mark Keppel High, 
Alhambra; Alan Johnston, Van Nuys High, Van Nuys; Harry Wat- 
kins, Inglewood High, Inglewood; Bennie Ford, Venice High, ‘Venice; 
Robert Lawrence, Wilson High, Los Angeles; David Johnston, 
Mark Keppel High, Alhambra; Lowell Doherty, Herbert Hoover 
High, San Diego; Bob Beavers, El Monte High, El Monte; Melvin 
Katz, Beverly Hills High, Beverly Hills; Phil Orville, Santa Monica 
High, Santa Monica; Gerald Parker, Venice High, Venice; Milo 
Webber, Pomona High, Pomona; and George Cooke, Hollywood 
High, Hollywood. 


Grades on the test ranged from a high of 85.4 to a low 
of 11.7, with a satisfactory distribution between these 
limits. A particularly good spread resulting in a strong 
“cream separator” action was obtained at the top of the 
list. 

In addition to the certificates of award, seven one-year 
college scholarships were available to the winners from 
the following schools: 
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California Institute of Technology, University of California (any 
campus), Occidental College, Pomona College, University of Red- 
lands, University of Southern California, and Whittier College. 


The winning students selected their scholarships from this 
list on the basis of rank and choice. 

The test this year was the best of those given so far, 
and represented a great deal of work on the part of the 
high school test committee. The committee consisted of: 


Professor Vernon L. Bollman, Occidental College, Chairman, Rev. 
Charles R. Coony, Loyola College, Professor Willard Geer, U.S.C., 
Professor Robert E. Holzer, U.C.L.A., Mr. Roy W. McHenry, 
Santa Monica City College, and Professor Kenneth Hurd, U.C.L.A. 
(member from host institution). 


The annual competitive physics test has become one of 
the Section’s most successful activities. The number of 
students and high schools participating has shown a very 
satisfactory increase from year to year, so that by now 
this annual test is regarded by both high school. physics 
teachers and their students as one of the important events 
of the year. The support of the colleges of this region, as 
evidenced by the scholarships offered in behalf of the 
Association, indicates the importance with which these 
institutions regard the test. 

FosTER STRONG, 
Secretary-Treasurer 





Oregon Section 


The eighteenth annual meeting of the Oregon Section 
of the American Association of Physics Teachers was held 
at the University of Washington on May 22, 1948. Seventy- 
two members and guests were present. Raymond T. 
Ellickson, President of the Section, presided. Fourteen 
papers were presented at the meeting. 

At the business meeting the following officers were 
elected by unanimous vote: President, Walter P. Dyke, 
Linnfield College; Secretary-Treasurer, Fred W. Decker, 
Oregon State College; Historian, Brother Godfrey Vassallo, 
University of Portland. The fall meeting will be held at 
Vanport College on November 13, 1948. 

1. High-vacuum technique—an undergraduate study. 
Harvey E. WEGNER, College of Puget Sound.—A high- 
vacuum system using a mercury diffusion pump and a 
gas injection system for ‘argon, helium, and neon has 
been constructed as an undergraduate project during the 
senior year. A few of the variety of tubes constructed for 
classroom demonstration are described, and some possible 
demonstrations for beginning physics students are dis- 
cussed. 

2. Philatelic physics. ALFRED B. BuTLER, Washington 
State College—At least 26 stamps have been issued bear- 
ing the portraits of physicists, and many more celebrate 
physical discoveries or applications. If the stamps are 
mounted individually with accompanying biographical 
sketch or description and are displayed in a hall showcase, 
a surprising number of students will be attracted and will 
linger to read the accompanying material. 

3. Factors determining teaching effectiveness in the 
elementary laboratory. T. W. LAsHor AND M. E. HoEHNeE, 


. 
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Reed College.—A preliminary report was given on a project 
for the compilation of a compendium of laboratory exer- 
cises. The report covered the first step: the establishment 
of a set of questions for the evaluation of exercises to be 
included in the compendium. 

4. The setting up of the mirror and lens equations and 
their unambiguous interpretation. A. E. HENNINGs, Uni- 
versity of British Columbia.—The wave front method is 
used to show the relations between , g, and f for mirror 
or lens. Real, inverted images are formed only when the 
wave front leaving the mirror or lens is converging. All 
other images are erect and virtual. Curves showing the 
relations between the object and image positions may be 
readily plotted. 

5. Microwave spectroscopy. A. VAN DER ZIEL, Uni- 
versity of British Columbia.—A survey of the techniques 
which are used in the investigation of rotational absorption 
spectra in the microwave region. The frequency of a klys- 
tron is swept linearly over about 10 Mc, and the signal 
generated by the klystron is passed through a wave-guide 
absorption cell and detected by a crystal diode. The de- 
tected signal is amplified and applied to the vertical de- 
flection plates of a cathode-ray oscillograph. On the hori- 
zontal plates a sweep voltage is applied which is syn- 
chronized with the frequency sweep of the klystron. In 
this way absorption spectra are displayed on the oscillo- 
graph screen. A microwave spectroscope built along the 
above principles is nearing completion at the University of 
British Columbia. 

6. Infra-red absorption spectra. A. M. Crookes, Uni- 
versity of British Columbia. 

7. Inertia effects in infra-red phosphors. RAYMOND T. 
Etickson, Reed College.—The visible light from an infra- 
red sensitive phosphor attains its maximum value after a 
period of the order of 0.001 second subsequent to the 
irradiation. An apparatus has been set up to measure this 
time of buildup. For one group of phosphors the time is 
independent of the intensity of the infra-red. For another 
group the build-up time is over one minute for very weak 
infra-red. 

8. Beta-ray energy measurements and disintegration 
schemes. K. C. MANN, University of British Columbia.— 
A discussion of some of the methods in the field of beta-ray 
spectroscopy to determine the energies of beta- and gamma- 
radiations from radioactive nuclei is given. This is ex- 
tended to the problem of determining the order of appear- 
ance of these radiations in the decay scheme. Some experi- 
mental results are quoted on Ra and Sb"™. 

9. Beta-decay at low energies. EUGENE P. Cooper, 
University of Oregon.—Indications of the failure of the 
Fermi theory of beta-decay are found in the very “rapid” 
decay of Hi’, the ‘‘abnormal” linearity with kinetic energy 
of the Cu® electrons at energies below 100 kev, and the 
apparent superabundance of Cu® positrons and Cd!” 
positrons in this low-energy region. It is suggested that 
radiation in the form of a continuous x-ray spectrum may 
accompany beta-decay in such a way as to be relatively 
most important at low energies where it could effectively 
modify the shape of the beta-spectrum. Such radiation 
may arise by virtue of “bremsstrahlung,”’ since the elec- 
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ANNOUNCEMENTS 


tron is decelerated in the field of the nucleus which emits 
it. A semiclassical estimate of this effect is given, and a 
method for obtaining a quantum-mechanical estimate is 
outlined. 

10. The humanities in the physics curriculum. RICHARD 
D. Murpay, c.s.c., University of Portland.—The trend to 
general education and the possible doubling of college 
enrollment by federal aid demand a broadening of the 
outlook of the physics teacher and the physics student. 
The need for the required integration will have to be met 
by a philosophy of life, which, in our present state, will 
have to be culled by the individual from the non-technical 
part of the curriculum. The faculty will have to bear the 
burden of fostering interests in the social problems of the 
day both by its counseling and its attitudes. 

11. Humanities and other nontechnical subjects in 
technical curricula. W. WENIGER, Oregon State College.— 
Attention is called to widespread activity in revising 
curricula with the purpose of offering students a more 
liberal education as well as preparing them for opportuni- 
ties in the newer fields of physics. The principal aim is to 
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keep alive an interest in all things, so that the physicists’ 
methods of attack may be applied in other fields. 

12, Cerenkov radiation in the vicinity of the threshold. 
Joun M. HarpinG and J. E. HENDERSON, University of 
Washington. 

13. Recent developments in the theory of spark dis- 
charge. RONALD GEBALLE, University of Washington. 

14, Scintillation counters. PHitip A. GOLDBERG, Unti- 
versity of Oregon.—A variety of problems associated with 
scintillation counter design, including such matters as 
phosphor selection, phosphor screen construction, photo- 
multiplier spectral response, signal-to-noise ratio, and a 
special light shield for the photo-multiplier tube are 
treated. Various optical accessories have been used be- 
tween the phosphor screen and the photo-multiplier tube 
to improve the energy differentiation characteristics of the 
counter. To increase the sensitive area of the counter, 
mirrors and a bank of photo-multiplier tubes have been 
used. 

WILLIAM R. VARNER, 
Secretary 


ANNOUNCEMENTS AND NEWS 


Meteorological Station at Brookhaven 
National Laboratory 


Two meteorological towers for weather observations in 
connection with nuclear research are being completed at 
Brookhaven National Laboratory. The taller tower, to be 
420 ft high, comparable in height to a 35- or 40-story build- 
ing, will be the tallest structure on Long Island. The 
second tower, already complete, is 160 ft high. At five 
levels on the small tower, and at eight levels on the tall 
tower, will be platforms where observers may take readings 
from meteorological instruments. 

Plans and design of the towers were supervised by the 
Brookhaven Meteorology Group, headed by Norman R. 
Beers, formerly Associate Professor of Aerological En- 
gineering at the Naval Post Graduate School, Annapolis, 
Maryland. Working with the group is a special station of 
the United States Weather Bureau, headed by Mr. 
Raymond C. Wanta. 

One of the chief purposes of the meteorological towers 
is to furnish information on wind currents in connection 
with safe operation of the Brookhaven nuclear reactor to 
be finished this fall. The nuclear reactor will be cooled by 
air. Large fans will conduct the air away from the pile, 
through an air duct, and up a 300-ft air stack on a hill 
near the reactor. In the air will be minute quantities of 
radioactive argon, a chemically inert gas, which will be 
dissipated harmlessly into the upper air. Study of wind 
velocity and direction, and atmospheric pressure, during 
operation will permit control of the pile so that air emanat- 
ing from the stack will not settle on or near the ground. 

Instruments on the weather tower and elsewhere on and 
near the Laboratory site will give precise information on 
wind currents. Operation of the reactor will be controlled 


at a low level or stopped entirely if it appears that radia- 
tion from the cooling stack under unfavorable conditions 
would not be completely harmless. With complete informa- 
tion on wind currents at upper and lower levels, and a 
knowledge of weather conditions on Long Island, scientists 
are thus assured that operation of the reactor will be 
harmless to neighboring communities, and to plant and 
marine life in the vicinity. 

Another feature of the tall weather tower will be a 
smoke stack to carry smoke specially created for the pur- 
pose of studying wind and weather. This is a 20-in. steel 
pipe designed to dissipate smoke produced by an army 
surplus M-1 smoke generator, similar to those used by 
naval vessels and beachhead troops in laying down smoke 
screens to conceal an operation. It produces a harmless oil 
fog smoke, 10,000 ft? per min, and at some time during 
tests it may create near the Laboratory a small cloud 
visible for several miles. 

Weather instruments on the towers will be mounted on 
beams that can be swung away from the towers, and then 
pulled in for observations. In addition, the observations 
will be recorded electrically on an instrument panel in a 
building 900 ft away. The building has to be a considerable 
distance from the towers so that eddies created by winds 
passing over the building will not disturb observations 
made by instruments on the towers. Two electrified cables 
between the towers will also carry instruments to record 
temperatures and wind differences between them. 

The meteorology group at the Laboratory will afford 
facilities for study and research by scientists of the per- 
manent staff and. by others who come on leave from their 
regular jobs elsewhere. The Weather Bureau will also be 
orienting and training U. S. Weather Bureau personnel, 
who will later be assigned to similar work elsewhere. 
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The following persons have been made members or 
junior members (J) of the American Association of Physics 
Teachers since the publication, Am. J. Physics 16, 369 
(1948), of the preceding list. 


Alishouse, Charles Clyde, Jr. (J), R. D. No. 1, Avonmore, Pa. 

Alpert, Nelson L., Rutgers University, New Brunswick, N. J. 

Anderson, Melvin Francis (J), Todd Union, University of Rochester, 
Rochester 3, N. Y. 

Baldwin, Leonard Bush, Jr. (J), 809 ‘“‘A’’ Birch Road, East Lansing, 
Mich. 

Banker, Eddie Keith (J), 2021 Home Ave., Columbus, Ind. 

Barlow, John Sutton, University of North Carolina, Chapel Hill, N. C. 

Beckerley, James G., P.O. Box 711, Mahwah, N. J. 

Belden, Ted F. (J), P. O. Box 5, Idyllwild, Calif. 

Bendler, Harry M., Michigan State College, East Lansing, Mich. 

Bielka, Robert Paul (J), 2211 No. 59th St., Seattle 3, Wash. 

Binkley, William Owen (J), 444 Summit Avenue, Hagerstown, Md. 

Blumenthal, Ralph H., 2939 Ocean Avenue, Brooklyn, N. Y. 

Bonner, T. W., Rice Institute, Houston, Tex. 

Bronson, S. Davis, 233 Wellington Rd. South, Garden City, Long 
Island, N. Y. 

Bryant, George Macon (J), 2025 Leighton Ave., Anniston, Ala. 

Bushkovitch, Alexander V., St. Louis University, 221 N. Grand Ave., 
St. Louis, Mo. 

Carter, David Southard (J), 6305 McCleery St., Vancouver, British 
Columbia, Can. 

Chu, En-Lung, Box 803, Stanford University, Calif. 

Clark, Donald Barnard (J), Box 876, Villanova, Pa. 

Cloyd, David Emerson (J), 6314 Manastash Place, Seattle 8, Wash. 

Clymer, Arthur Benjamin, MIT Graduate House, Cambridge, Mass. 

Crutchfield, Bobby L., University of Texas, Austin, Texas. 

Dahl, Adrian H., 114 Disston Road, Oak Ridge, Tenn. 

Dam, Cecil F. (J), 450 Egleston Avenue, Kalamazoo 30, Mich. 

Davidson, Kenneth W. (J), Rt. No. 2 Abernathy, Tex. 

Day, Franklin, Box 579, Pikeville, Ky. 

Deake, Lt. Standish, Woodstock Country School, Woodstock, Vt. 

DePangher, John De, 2103 Grant Street, Berkeley 3, Calif. 

Di Bella, Angelo (J), 447 North Alameda, Las Cruces, N. M. 

Dick, Louis A., Jr., (J), 627 South Hastings Ave., Hastings, Nebr. 

Doller, Dr. Samuel, 1700 Crotona Park East, Bronx 60, N. Y. 

DuBois, James G. E., 11272 Nebraska Ave., West Los Angeles 25, Calif. 

Everhart, Edgar, 1-B Gibson Terrace, Cambridge, Mass. 

Fairweather, Stephen H., 1088 Fairfield Ave., Bridgeport, Conn. 

Falkoff, David Leo, University of Notre Dame, South Bend, Ind. 

Fankuchen, Isidor, Polytechnic Institute of Brooklyn, Brooklyn 2, N. Y. 

Ferguson, George Alonzo, 5120 Just Street, N. E., Washington 19, D. C. 

Fink, Dr. George A., 147-30 Barclay Ave., Flushing, Long Island, N. Y. 

Fuguson, Mary C., University of Nevada, Reno, Nev. 

Fuson, Nelson, Johns Hopkins University, Baltimore 18, Md. 

Gabriel, Walter James (J), 2037 Windsor Place, Fort Worth 4, Tex. 

Galginaitis, Simeon V., University of Louisville, Louisville, Ky. 

Gerjuoy, E., University of Southern California, Los Angeles, Calif. 

Gogarn, Frances E., 634 S. State St., Big Rapids, Mich. 

Goodrich, Raymond F. (J), 6500 Newgard Ave., Chicago 26, Ill. 

Goodson, Dale (J), William Jewell College, Liberty, Mo. 

Grafa, Julian B., 121 S. Knoblock, Stillwater, Okla. 

Gray, Donald, University of Southern California, University Park, 
Los Angeles 7, Calif. 

Greenwood, Joseph P., 2539 Bonnie Lane, Stockton, Calif. 

Harlor, T. Lyle (J), 176 East Burlington Road, Riverside, Ill. 

Harrison, Charles Edward, 4530a West Poplin St., St. Louis 10, Mo. 

Harrison, Melvin A. (J), University of Washington, Seattle 5, Wash. 

Harrison, William Hartley, 6363 Santa Monica Blvd., Hollywood 38, 
Calif. 

Hart, Hiram, 510 W. 188 Street, New York 33, N. Y. 

Hatke, Sister M. Agnes, St. Francis College, 2701 Spring St., Fort 
Wayne, 8, Ind. 

Heinrich, Edwin P., 6204 Macarthur Blvd., Washington 16, D. C. 

Heirtzler, James R., Box 6941, University Station, Baton Rouge 3, 
La. 
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New Members of the Association 


Herman, Eugene B. (J), 120 East Foster Ave., State College, Pa. 

Holter, Norman J., 113 North Main Street, Helena, Mont. 

Hoyle, Hughes Bayne, Jr., 1907 Vail Ave., Apt. 3, Charlotte, N. C. 

Hsiang, Jen-Sen, Rutgers University, New Brunswick, N. J. 

Hueston, Roland A., Jr., 11 Lindberg Road, Framingham, Mass. 

Hughes, Darrell Stephen, 801 Park Blvd., Austin, Tex. 

Humphrey, Samuel S. (J), Box 146, Storrs, Conn. 

Hurlbut, James Franklin (J), 1401 Fairfax Street, Apt. 1, Denver, Colo. 

Johnson, John Lowell (J), Montana State College, Bozeman, Mont. 

Johnson, Walter E., Benedict College, Columbia, S. C. 

Kennel, John Maurice (J), 304 N. Bishop St., Oxford, Ohio. 

Koda, Nobuo John (J), 1919 W. 44th St., Cleveland, Ohio. 

Lance, Harvey W., Apt. 308, 201 Elmira St. S. W., Washington 20, 
B.C. 

Landis, Harry M. (J), 1909 Northwood Blvd., Royal Oak, Mich. 

Lenahan, John James (J), 605 Lafayette Ave., Buffalo 9, N. Y. 

Lindsey, Tate J., Box 172, Clemson, S. C. 

MacCracken, Elliott B., 262 Hargadine St., Ashland, Oregon. 

Masters, Joseph (J), 20 West Street, Hanover, N. H. 

Maxwell, Arthur E. (J), Box 383, State College, N. M. 

McFarland, Charles E. (J), 1718 Janet Place, Kirkwood 22, Mo, 

McMichael, Harvey G. (J), 1007 E. Marion St., Seattle 22, Wash. 

Mehl, Clarence Robert (J), 2710 Third Avenue North, Great Falls, 
Mont. 

Miller, James Richard, 909 Prospect Ave., Takoma Park, Md. 

Mitchell, Robert C. (J), Box 411, State College, N. M. 

Morris, Thomas P. (J), R. F. D. 1, Liberty, Missouri. 

Naiditch, Sam, Duquesne Univ., Pittsburgh 19, Pa. 

Noyes, Bradford, 4087 N. Illinois St., Indianapolis, Ind. 

Oerlein, Karl F., Army Chemical Center, Chemical Corps School, Md. 

Ohralik, Karl A. (J), 2182 West 105 St., Cleveland, Ohio. 

Page, Thornton Leigh, University of Chicago, Chicago 37, Ill. 

Parke, Nathan Grier III, Spencer Brook Road, Concord, Mass. 

Pepper, Robert Richard, Purdue University, W. Lafayette, Ind. 

Puechl, Karl H. (J), S. Londonderry, Vt. 

Rainey, Richard Owen (J), 36 Veterans Village, University of Red- 
lands, Redlands, Calif. 

Rasco, William Allen, University of Texas School of Dentistry, Hous- 
ton, Tex. 

Rawson, Edward B. (J), Route 5, Frederick, Md. 

Reed, Benjamin L. III, 4402 La Cara, Long Beach 4, Calif. 

Rider, John Allen, Jr., Box 276, Southwestern Louisiana Institute, 
Lafayette, La. 

Rieu, Emile Henry (J), 1018 Kirkpatrick Street, North Braddock, Pa. 

Rosario, Dr. Leticia Del, University of Puerto Rico, Rio Piedras, 
Puerto Rico. 

Ruble, George B. (J), 28 Converse Ct., Burlington, Vt. 

Schiff, Daniel (J), 186 S. Cassady, Columbus 9, Ohio. 

Scibor-Marchocki, Romuald Ireneus, General Delivery, Pasadena, 
Calif, 

Seely, Horace E., Jr. (J), 941 E. 83 St., Chicago 19, Ill. 

Singleton, Burt Newman, Jr. (J), 1404 W. Palmetto St., Florence, S. C. 

Slack, Lewis, Washington University, St. Louis 5, Mo. 

Stanley, Charles Vincent (J), 124 Buckeye Street, Hamilton, Ohio. 

Stolov, Harold L., 100 VanCortlandt Park South, Bronx 63, N. Y. 

Struwin, Richard, 284 Champion St., Battle Creek, Mich. 

Thielman, Charles E. (J), 303 12th St., S.E., Paris, Texas 

Thompson, Ross D. F., Rockefeller Hall, Cornell University, Ithaca, 
N.Y. 

Van Bergen, Francis E., 15 Spoorweglaan, St. Niklaas-Waas (Belgium). 

Walker, John Pearson, Jr., 913 W. College, State College, Pa. 

Washington, Mauree L., 1825 Vernon St., N. W., Washington 9, D. C. 

Wassink, H. J., R. No. 5, Holland, Mich. 

Webb, Herbert Godfrey (J), 1002 S. Davis, Kirksville, Mo. 

Willardson, Robert K. (J), 329 Callister Court, Provo, Utah. 

Willett, Donald Kenneth (J), R. F. D. 2, Concord, Mass. 

Winter, Rolf G. (J), General Delivery, Carnegie Inst. of Tech., Pitts- 
burgh 13, Pa. 

Woodruff, Willis P., Jr. (J), 306 Woodland Ave., Birmingham 9, Ala. 

Zechiel, Leon Norris (J), 2245 Neil Ave., Columbus 1, Ohio. 
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